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THE  SIMULATION  LANGUAGE  CSMP-IO(ARL) 


N.  E./gILBERT  ami  P.  G,  NANKIVELL 


SUMMARY 

A description  is  given  of  the  'block  oriented'  simulation  language  CSMP-!0(ARL). 
which  has  been  developed  from  CSMP-IO  and  is  written  mainly  in  FORTRAN  IV  for  a 
P DP-10  computer.  The  major  improvement  made  has  been  to  incorporate  'user-defined' 
blocks,  which  are  written  as  FORTRAN  subroutines.  A large  number  of  outputs  may  be 
defined  within  these  subroutines  by  using  'dummy'  or  'user  output'  blocks.  Two  other 
major  improvements  have  been  made,  both  enabling  the  .saving  of  appreciable  core  storage. 
Firstly,  arrays  necessary  to  store  information  on  each  block  u.sedarc  automatically  expanded 
to  the  size  determined  by  the  u.ser.  Secondly,  the  language  is  divided  into  a modelling 
program,  which  is  used  to  perform  the  model  simulation  and  store  the  output  in  a binary 
file  on  a specified  storage  device,  and  an  output  program,  which  is  used  to  print  and  plot 
the  character  conversion  of  the  binary  file.  ^ 
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1.  INTRODUCTION 

Programming  of  both  continuous  and  discrete  dynamic  systems  for  digital  computers  has 
been  aided  considerably  in  recent  years  by  the  development  of  simulation  languages.  These 
languages  enable  the  user  to  make  a relatively  simple  transition  from  the  environment  of  an 
analogue  computer  to  that  of  a digital  computer.  Simulation  languages  have  advantages  over 
high  level  programming  languages  such  as  FORTRAN  and  ALGOL  of  economy  in  program- 
ming time  through  the  sorting  of  statements  and  ease  of  modification,  use  of  centralised 
integration,  and  specialised  output  facilities.  There  are  two  types  of  simulation  language,  ‘block 
oriented'  and  ‘equation  oriented’;  the  latter  resemble  a high  level  programming  language  such 
as  FORTRAN,  while  the  former  are  e.xpressed  in  coded  form  with  the  aid  of  a block  diagram.* 
The  interactively  controlled  ‘block  oriented'  language  CSMP-10'  became  available  for  use  on 
the  POP-10  computer  at  ARL  but,  because  of  its  inability  to  handle  complex  algebraic  expres- 
sions, was  of  limited  usefulness.  Consequently,  a modified  version,  called  CSMP-IO(ARL),  has 
been  developed  whose  major  improvement  is  the  incorporation  of  user-defined  blocks  which 
are  written  as  FORTRAN  subroutines.  CSMP-10  was  developed  from  a series  of  languages, 
namely  (in  order  of  deselopment)  PACTOLUS^  for  an  IBM-1620,  1130-CSMP’'^  for  an  IBM- 
1 130,  and  CSMP-9^  for  a PDP-9.  Besides  CSMP-I0(ARL),  there  have  been  a number  of  other 
off-shoots  with  extended  capabilities,  including  ll.fOCSMP**  ’ written  for  an  IBM-1130,  and 
allowing  the  use  of  five  short  user-defined  FORTRAN  subroutines  and  several  other  improve- 
ments. Another  version,  ***CSMP,"  offers  several  improvements  which  also  appear  in 
CSMP-I0(ARL),  but  allows  only  two  short  user-defined  subroutines  written  in  BASIC. 

An  outline  description  of  CSMP-I0(ARL)  is  given  in  Section  2,  which  includes  also  the 
improvements  made  to  CSMP-10.  Section  3 shows  how  a specific  mathematical  model  may  be 
represented  using  various  statements,  while  Section  4 gives  the  Teletype  responses  necessary  for 
the  beginner  to  complete  the  simulation  of  his  specified  mathematical  model.  Details  of  more 
advanced  features  of  the  language  are  given  in  Section  5 to  9.  In  addition  to  the  many  changes 
made  to  CSMP-10,  the  Teletype  messages  and  corresponding  responses  have  been  completely 
revised,  and  much  greater  detail  is  provided  here  than  in  descriptions  of  previous  versions  of 
the  language. 

The  language  CSMP-IO(ARL)  is  currently  being  used  at  ARL  in  the  computer  simulation 
of  the  behaviour  of  helicopters  while  under  manual  or  automatic  flight  control. 


• A block  diagram  (e  g.  Fig.  2)  consists  of  a number  of  linked  modules  (called  blocks),  each  one 
representing  a particular  function  or  operation. 
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2.  OUTLINE  DESC  RimON  OF  CSMP-IO(ARL) 

2.1  General 

CSMP-10^ARL)  is  written  mainly  in  FORTRAN  IV  for  the  PDP-IO  FORTRAN  compiler 
‘F40  \'ersion-27‘  running  under  the  operating  system  ‘FOROTS’.  The  language  consists  of  a 
modelling  program  named  BOMMP  (Block  Oriented  Mathematical  Modelling  Program)  and 
an  output  program  named  TRANS  (Translation).  In  conjunction  with  various  statements, 
control  parameters,  and  FORTRAN  subroutines  defined  by  the  user  in  order  to  represent  a 
specific  mathematical  model,  the  modelling  program  is  used  to  perform  the  simulation  and 
store  the  output  in  a file  in  binary  form  on  a specified  storage  device.  The  output  program  is 
then  used  to  print  and  plot  the  character  conversion  of  the  binary  file  in  either  tabular  or  graphical 
form.  The  output  program  may  be  used  independently  of  the  modelling  program  as  a general 
purpose  output  program  (provided  data  are  supplied  in  the  appropriate  format).  Only  a brief 
description  of  its  use  is  therefore  included  here  (in  Section  4.2);  a complete  description  is  given 
in  Reference  9.  Although  the  modelling  program  is  designed  to  be  run  interactively  from  a 
Teletype,  it  may  also  be  run  non-interactively  during  Batch  or  On-Line  processing.  In  the 
modelling  program,  the  subroutines  CPU,  EXPAND,  FCFIECK,  and  TTYCHK  are  written  in 
MACRO-IO.  the  assembly  language  for  the  PDP-IO. 

The  modelling  program  allows  a specific  mathematical  model  of  a dynamic  system  to  be 
letiresented  by  control  parameters  and  three  types  of  statements,  viz,  configuration,  parameter, 
and  function  statements.  The  configuration  statements  describe  the  block  structure  of  a particular 
model:  each  statement  consists  of  a block  number,  a block  type  (e.g.  integrator,  adder),  and  a 
list  of  which  other  blocks  (up  to  three)  supply  the  block  with  input  signals.  The  parameter  state- 
ments specify  numerical  values  of  parameters  associated  with  the  configuration  statements, 
while  function  statements  specify  coordinate  pairs  used  to  generate  a function.  The  configuration 
statements  may  be  listed  in  any  order;  the  modelling  program  (BOMMP)  sorts  them  into  an 
executable  order.  Following  these  statements,  integration  and  output  control  parameters  are 
specified.  Integration  is  performed  numerically  using  a second  order  Runge  Kutta  (i.e.  Modified 
Euler)  method  with  a fixed  integration  interval.  The  user  specifies  the  lower  and  upper  time 
limits  and  the  interval  thus  giving  equispaced  time  steps  at  which  the  calculations  are  performed 
and  required  output  values  are  stored.  For  the  integration  method  used  (see  Section  6.2),  there 
are  two  stages  of  calculation.  The  second  stage  provides  a solution  of  the  simulation  at  the  present 
time  step,  while  the  first  stage  provides  a solution  at  a time  value  midway  between  the  present 
and  previous  time  steps.  This  'mid-point'  time  value  is  referred  to  as  a 'half  time  value'  or 
'half  time  step'. 

2.2  Improvements  to  C'S\1P-I0 

(o)  ( • scr-definvd  block  .s 

Because  of  its  in;ibility  to  handle  complex  algebraic  expressions  cither  directly  in  the 
configuration  statements  or  indirectly  in  user-defined  subroutines  written  in  FORTRAN, 
(SMP-IO  cannot  be  considered  a very  useful  general  purpose  simulation  language.  The  major 
improvement  made  in  CSMP-IO(ARL)  therefore  has  been  to  incorporate  user-defined  blocks, 
which  arc  written  as  F ORTRAN  subroutines  and  are  compiled  and  loaded  into  core  with  the 
compiled  version  of  BOMMP.  The  current  version  of  BOMMP  allows  a maximum  of  fifteen  of 
these  blocks;  however,  this  number  may  be  increased  indefinitely  within  the  limits  of  the  com- 
puter core  storage  by  simple  coding  alterations  to  BOMMP.  These  blocks  may  have  only  three 
sorted  block  inputs  (as  declared  in  the  configuration  statement),  but  may  have  any  number  of 
unsoricd  block  inputs.  'Dummy'  or  ‘user  output  olocks  have  been  created  whose  output 
value  is  defined  in  the  appropriate  user-defined  subroutine,  thus  allowing  a large  number  of 
additional  outputs  for  each  of  these  blocks. 

(6)  l.'xpoiiddhlc  arrow 

With  CSMP-10.  the  total  number  of  blocks  is  limited  to  seventy-five,  which  includes  a 
maximum  of  twenty-five  integrator  and  unit  delay  blocks,  and  three  function  blocks.  These 
limits  could  be  increased  by  resetting  the  appropriate  DIMENSION  statements  in  BOMMP. 
but  to  do  so  would  reserve  a large  amount  of  core  storage  that  may  not  be  required  for  most 
mi’dels  Therefore,  the  facility  has  been  introduced  in  CSMP-IO(ARL)  whereby  the  user  specifies 
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(i)  the  maximum  block  number, 

(ii)  the  total  number  of  integrator  and  first  order  lag  blocks, 

(iii)  the  number  of  unit  delay  blocks,  and 

(iv)  the  number  of  function  blocks. 

The  appropriate  arrays  in  BOMMP  are  then  automatically  expanded  to  the  necessary  size  so 
that  a minimum  of  core  storage  is  used.  The  provision  of  this  facility  required  the  inclusion  of  a 
subroutine  named  EXPAND,  which  as  previously  stated  is  written  in  the  assembly  language 
MACRO- 10  for  the  PDP-10.  For  other  computers,  this  subroutine  would  have  to  be  rewritten  in 
the  appropriate  language;  alternatively,  the  facility  may  be  removed  and  array  limits  set  as  in 
CSMP-IO. 

(<■)  Block  output 

In  CSMP-IO,  a significant  amount  of  core  storage  is  used  to  code  the  output  section  of  the 
program.  By  dividing  the  language  into  a modelling  program  and  an  output  program,  the  core 
storage  required  for  coding  is  significantly  reduced.  The  size  of  the  problem  that  can  be  handled 
in  CSMP-IO(ARL)  is  then  determined  by  the  requirements  of  the  modelling  program  alone.  A 
large  number  of  improvements  to  the  output  part  of  the  language  have  been  effected  without 
reducing  the  size  of  the  problem  that  can  be  handled.  Also,  because  the  'raw'  output  data  are 
stored  in  a file  in  binary  form  by  the  modelling  program,  various  forms  of  output  of  the  same 
data  may  be  obtained  without  repeating  the  simulation.  As  a further  saving,  an  output  value  is 
only  stored  if  it  changes  by  more  than  a specified  percentage  from  the  last  value  stored  at  a previous 
time  step. 

(d)  Other  Improxements  to  CSMP-IO 

In  addition  to  the  three  major  improvements  described  above,  a large  number  of  minor 
improvements  have  been  made,  many  of  them  being  corrections  to  errors  in  the  available  version 
of  CSMP-IO.  The  significant  improvements  include: 

(i)  an  improved  iterative  method  of  solving  an  implicit  equation; 

(ii)  improved  integration  facilities  which  include  a first  order  lag  block  and  ‘hold’  and 
‘reset’  operations  for  the  integrator  block; 

(iii)  a function  block  with  any  number  of  arbitrarily  spaced  abscissae  and  where,  in  addition 
to  linear  interpolation,  either  a polynomial  of  degree  two  or  higher,  or  a user-defined 
function,  may  be  used  to  calculate  the  function  value; 

(iv)  a debugging  facility; 

(v)  alphanumeric  labelling  of  each  block;  and 

(vi)  the  capability  of  continuing  execution  of  the  model  at  a later  stage. 

3.  MODEL  REPRESENTATION 

To  represent  a dynamic  system  in  CSMP-I0( ARL),  the  user  should  first  produce  a block 
diagram  of  the  type  suitable  for  programming  an  analogue  computer.  The  list  of  available 
block  functions  in  CSMP-I0(ARL)  is  given  in  Table  1 (see  Section  3.1  for  associated  notation). 
More  complex  mathematical  operations  may  be  defined  by  the  user  as  user-defined  blocks 
written  as  FORTRAN  subroutines.  The  block  diagram  should  then  be  translated  into  configura- 
tion, parameter,  and  function  statements.  These  statements,  which  are  outlined  briefly  in  Section 
2.1.  arc  described  in  detail  below  in  Section  3.1.  To  complete  the  model  representation,  the  user 
needs  to  specify  integration  and  output  control  parameters  (see  Sections  2.1  and  4.1.6). 

In  order  to  understand  many  of  the  features  of  the  language,  particularly  the  more  advanced 
ones,  it  is  important  to  know  the  algorithm  for  sorting  the  configuration  statements  into  an 
executable  order  and  the  processes  by  which  the  model  is  executed.  These  two  processes  are 
therefore  described  below  in  Sections  3.2  and  3.3  respectively. 

3.1  Types  of  Stafements 

As  outlined  in  Section  2,1,  there  are  three  types  of  statements  associated  with  the  model 
representation,  viz,  configuration,  parameter,  and  function  statements.  The  notation  associated 
with  these  statements  and  the  integration  parameters  is  defined  as  follows  (see  Fig.  1 for  corres- 
ponding block  diagram  notation); 


B 

Bl,  B2,  B3 
M 

PI,  P2,  P3 
T 

X 

XI,  X2,  X3 
h 

t 

^0 


block  number  (i.e.  an  integer  identifier) 

blocks  which  supply  input  signals  to  a given  block  (referred  to  as  inputs) 
integer  parameter  used  only  with  the  block  types  F,  I,  T,  Tl,  V,  and  Y ; zero  for 
all  other  blocks  (see  Table  I and  appropriate  sections) 
parameters  associated  with  a given  block 

block  type  (library  given  in  Table  I ; each  type  has  an  associated  reference 
number) 

output  value  of  block  B 

output  value  of  blocks  Bl,  B2,  and  B3  respectively 

lime  interval  for  integration 

time 

initial  time 


Errors  and  their  diagnostic  messages  relating  to  the  above  specification  statements  are 
presented  in  Appendix  B. 


TABLE  1 
Block  Type  Library 

(refer  Appendix  A) 


Name 

Type 

Type  No, 

No,  of 
Inputs 

No  of 
Parameters 

Description 

Bang  bang 

B 

2 

1 

0 

X = 1 if  XI  >0 
= -1  ifXI  ■ 0 

Dead  space 

D 

4 

1 

2 

X = XI-PI  if  XI  > PI 
= XI -P2  if  XI  ■ P2 
= 0if  P2  < XI  « PI 

Function 

F 

6 

1 

1 

X = f(X  1 ).  Interpolation  from 
a function  generated  by  co- 
ordinates with  arbitrarily- 
spaced  abscissae;  linear  if 
PI  = 0 or  1,  of  degree  PI  if 
PI  > 1,  user-defined  if  PI  = 
— 1,  —2,  or  —3,  Uses  P3  and 
M internally  (see  Section  3. 1 .3) 

Gain 

G 

7 

1 

I 

X = P1*XI 

Flalf  power 

H 

8 

1 

0 

X = x'XI 

Integrator 

1 

9 

3 

3 

X = P1  . jljdt 

where  y = Xl(l  - P2)  • P3; 
if  X2  ^0;  RESET  TO  PI  if 
X3  ^ 0,  Uses  M internally 
(see  Section  6,2). 

Jitter 

J 

10 

0 

0 

X = {random  number  uni- 
formly distributed  between 
— 1 and  3 1} 

Constant 


K 


0 


X = PI 


TABLK  1 — continued 


Name 

Type 

Type  No. 

No.  of 
Inputs 

No.  of 
Parameters 

Description 

Limiter 

L 

12 

1 

2 

X PI  ifXl  > PI 
= P2  if  XI  P2 
XI  if  P2  ^ XI  $ PI 

Magnitude 

M 

12 

1 

0 

X X 

Negative 

clipper 

N 

14 

1 

0 

X XI  if  XI  > 0 
= 0 if  X 1 ^0 

Olfsct 

() 

15 

1 

1 

X = XI  PI 

Positive 

clipper 

P 

16 

1 

0 

X XI  if  XI  ■ 0 
0 if  XI  >0 

Quit 

Q 

17 

2 

0 

Terminates  run  if  XI  X2; 

types  message 

RL'X  TRKMI.WirED  BY  A 
Q BLOCK 

Relay 

R 

18 

,2 

0 

X X2  if  XI  > 0 
X3  if  XI  0 

Time  pulse 
generator 

T 

20 

1 

1 

Geni  ates  pulse  train  of  period 
PI  starting  when  XI  ^ 0. 

Lises  P2  and  M inlernally  (see 
Appendix  A;  also.  Section  8.2 
for  use  in  discrete  system 
simulation) 


f irst  order 

T1 

.20 

2 

1 

Solution  of  the  first  order 

lag 

dilferential  equation 

,\  P2,/\./r  XI  -X2  X.1 
vs  here  X PI  when  t r„. 
Uses  M internally  (see  Section 


6,.'<) 

I' nit  delay  U 21  I I X PI  when  i r„ 

X I at  r h 2 when  i ■ la- 
Uses  P2  internally  when  the 
input  Bl  is  a U block  (see 
Section  6.1;  also  Section  8.2 
for  use  in  discrete  system 
simulation) 

UO  .1  I 0 X is  defined  in  the  user- 

defined  subroutine  associated 
with  the  input  Bl  (see  Section 
7) 


User  output 


TABLE  1 — continued 


Name 

Type 

Type  No. 

No.  of 
Inputs 

No.  of 
Parameters 

De.scription 

User- 

defined 

Urn 

30^  m 

3 

3 

X is  defined  by  the  user  in 
subroutine  USERm,  where  m 
is  an  integer,  1 to  15  (see 
Section  7) 

Vacuous 

V 

22 

0 

1 

Use  with  wye  block  to  solve 
an  implicit  equation  (e.g. 
y = /(y)).  PI  is  the  initial 
approximation  of  y,  i.e.  y„. 
Uses  P2,  P3,  and  M intern- 
ally (see  Section  8.1) 

Weighted 

summer 

W 

23 

3 

3 

X = P1*XI  • P2*X2  • P3*X3 

Multiplier 

X 

24 

2 

0 

X = XI»X2 

Wye 

Y 

25 

2 

2 

Used  with  vacuous  block  to 

solve  an  implicit  equation 
(e.g.  V = /'(y)).  X = y:  PI  is 
the  error  tolerance  for  con- 
vergence, Ec;  P2  is  the  maxi- 
mum number  of  iterations, 
Ni;  if  not  set,  Ni  is  assumed 
to  be  20,  Uses  P3  and  M 
internally;  also,  dual  use  of 
P2  (see  Section  8.1) 


Zero  order 
hold 

Z 

26 

2 

1 

X = PI  when  t = r„ 

= XI  if  X2  > 0 when 
t '■  t„:  X is  unchanged  if 
X2  ^ 0 when  / > r„:  P2  is 
used  internally  to  store  X at 
the  previous  half  time  step 
(see  Section  8.2  for  use  in 
discrete  system  simulation) 

Summer 

27 

3 

0 

X = ±XI±X2±X3.  The 

inputs  Bl,  B2.  B3  may  be 
preceded  by  a plus  or  a minus 
sign  in  the  configuration  state- 
ment 

Divider 

/ 

29 

2 

0 

X = XI, 'X2 

Inverter 

— 

28 

1 

0 

X = -XI 

Blank 

5 

0 

0 

Causes  the  previous  state- 
ment with  the  same  block 
number  to  be  deleted  but  not 
replaced 

3.1.1  Configuralion  Statements 

The  eonliguralion  statements  describe  the  block  structure  of  a model  and  have  the  form 


B.  T.  Bl.  B2,  B3;  "label”  S "comment" 

The  arguments  arc  separated  by  commas  (as  shown)  or  blank  characters.  The  modelling 
program  reserves  block  number  1.  whose  output  value  is  the  current  time,  /.  A blank  block 
number  (i.c.  two  consecutive  ‘carnage-returns')  terminates  the  set  of  configuration  statements. 
The  library  of  available  block  types  is  given  in  Table  I.  The  following  cvampics  of  configuration 
statements  illustrate  the  various  ways  these  statements  may  be  typed  (evpianations  given  below  );t 

2 1)  16  r ALPHA  S IN  RADIANTS 

1,  4.  5,  6 SACCEI.KRATION 

7.  K;  PHI 

8 \V  12  •)  B 3 is  blank  so  that  X?  is  assumed  zero 

The  inputs  Bl.  B2.  and  B3  must  be  valid  block  numbers,  or  zero  or  blank  to  signify  the 
absence  of  an  input.  If  the  block  type  is  a ‘summer’  (i.e.  T is  ‘ • ‘).  the  inputs  Bl.  B2,  and  B3 
m.iy  be  negative  integers.  If,  for  c.xample.  Bl  equals  —3.  this  means  that  the  output  value  of 
block  number  3 will  be  multiplied  by  —I  when  used  as  the  input  signal  to  the  given  ‘summer’ 
block.  Although  the  general  delinition  of  a block  type  may  assume  the  presence  of  specific 
inputs,  one  or  more  of  the.se  inputs  may  be  left  blank,  in  which  case  the  corresponding  value  of 
each  of  these  inputs  is  assumed  to  be  zero  (e.g.  in  the  above  examples,  B3  for  block  number 
8).  Block  types  that  often  do  not  require  all  their  specifiable  inputs  are  ‘integrator"  (I),  ‘first 
order  lag’  (Tl),  ‘user-defined’  (Urn),  ‘weighted  summer'  (W),  and  ‘summer’  ( ). 

The  semi-colon  in  the  configuration  statement  is  optional:  when  used,  the  ten  characters 
following  may  be  used  to  provide  an  alphanumeric  label  for  B,  which  is  retained  when  storing 
the  model  (see  Section  4.1.5)  and  is  used  to  label  the  printed  and  plotted  output.  The  dollar  sign 
terminator  is  also  optional;  all  characters  following  it  are  ignored  and  may  therefore  be  used  as 
comments. 

The  configuration  statements  may  be  input  in  any  order.  The  sort  process,  described  in 
Section  3.2.  determines  the  order  in  which  the  statements  are  executed.  If  the  user  types  a con- 
figuration statement  having  a block  number.  B.  that  has  already  been  used  as  the  block  number 
of  a previous  statement,  the  previously  defined  statement  is  replaced  by  the  new  statement  unless 
the  type,  T,  of  the  new  statement  is  blank,  in  which  case  the  previous  statement  is  deleted  but  not 
replaced  (see  examples  in  Section  4,1.4). 

3.1.2  Parameter  Statements 

F’aramctcr  statements  specify  parameter  values,  e.g.  initial  conditions  and  multiplication 
constants,  associated  with  the  configuration  statements  (see  Table  I).  and  have  the  form 

B.  PI.  P2.  P3  S "comment" 

As  with  configuration  statements,  the  arguments  may  be  separated  by  eommas  or  blanks,  and 
two  consecutive  ‘carriage-returns'  terminate  the  set  of  parameter  statements.  Here  again,  the 
dollar  sign  followed  by  a comment  is  optional,  and  previous  parameter  statements  may  be 
replaced  by  repeating  the  block  number.  The  parameters  PI,  P2.  and  P3  arc  associated  with  the 


t Teletype  messages  typed  by  the  user  are  shown  in  bold  upper  case;  messages  typed  by  the 
appropriate  computer  program  arc  shown  in  italic  uppercase.  Comments  on  Teletype  messages 
are  shown  in  upper  and  lower  ease  print  alongside  the  appropriate  message. 
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block  number  B appearing  in  the  configuration  statements.  For  example,  when  the  parameter 
statement 


21,  1.5  S INITIAL  C ONDITION  FOR  \ 

is  typed  for  a model  with  the  previously  defined  configuration  statement 

21.  I.  lO.l 

the  initial  value  of  block  number  21  is  .set  equal  to  1.5  (as  defined  b>  the  specification  for  an 
integrator  block  type  in  Table  1). 

Aecause  storage  for  the  parameters  is  allocated  for  each  of  the  blocks  defined  in  the  configu- 
ration statements,  any  parameters  not  set  will  assume  a value  of  zero.  To  delete  the  parameters  of 
a particular  block,  the  user  sets  the  parameters  equal  to  zero  by  typing  only  the  block  number. 

3.1.3  Function  Statements 

Function  statements  specify  coordinate  pairs  used  to  generate  a function  for  each  block  of 
type  ‘F’  defined  in  the  configuration  statements.  The  appropriate  block  number  is  first  typed 
following  the  Teletype  message 

BLK  NO. 

and  each  coordinate  pair,  i.e.  abscissa  followed  by  ordinate,  is  then  typed  on  a separate  line 
following  the  message 

COORD  P.^IRS: 

(see  example  in  Section  4.1.4).  Two  consecutive  ‘carriage-returns’  terminate  the  list  of  co- 
ordinate pairs.  The  first  of  the  two  Teletype  messages  above  is  then  repeated  so  that  further 
functions  may  be  specified:  a blank  block  number  (i.e.  ‘carriage-return')  terminates  the  set  of 
function  statements. 

The  coordinate  pairs  define  a function  with  arbitrarily  spaced  abscissae.  Between  each 
adjacent  pair  of  abscissae,  the  value  of  the  function  is  obtained  by  interpolation.  The  inter- 
polation is  linear  if  PI  = Oor  I.  of  degree  PI  if  PI  ■ 0.  and  user-defined  if  PI  —1,  —2,  or  3 
(see  Appendix  C).  Outside  the  abscissae  limits,  the  function  value  is  obtained  by  extrapolation. 
The  coordinate  pairs  may  be  typed  in  any  order;  on  typing  each  pair,  the  appropriate  list  is 
rearranged  such  that  the  abscissae  are  in  order  of  increasing  magnitude.  \V  hile  entering  the  list, 
coordinate  pairs  may  be  deleted  or  replaced  (see  example  in  Section  4.1.4).  The  list  may  also 
be  modified  at  subsequent  stages  (see  Section  5.3).  The  internal  use  of  the  parameter  P3  and 
integer  parameter  M in  locating  the  position  of  the  data  for  each  function  (F)  block  in  the 
appropriate  expandable  array  (used  to  store  all  the  function  coordinate  pairs)  is  described  in 
Appendix  C. 

For  a polynomial  of  degree  PI.  at  least  PI  • 1 coordinate  pairs  must  be  defined.  There  is  no 
limit,  subject  to  sufficient  core  storage  being  available,  to  the  number  of  coordinate  pairs  that 
may  be  defined  for  each  function.  Flowever,  the  number,  n say,  of  F blocks  specified  when 
dimensioning  the  appropriate  expandable  array  (see  Section  4.1.2)  must  not  be  less  than  the 
actual  number  of  F blocks  used,  and  the  total  number  of  coordinate  pairs  of  all  F blocks  must 
not  exceed  1 5/i  (i.e.  to  exceed  an  average  of  1 5 points  per  F bliK'k.  n must  be  made  larger  than  the 
actual  number  of  F blocks  used). 

3.2  Sort  Process 

The  purpose  of  the  sort  process  is  to  arrange  the  blocks,  represented  by  their  configuration 
statements,  into  an  order  such  that  for  each  statement  the  input  signals  have  been  previously 
defined.  The  list  of  blocks  obtained  in  this  way  is  termed  the  ‘sort  list’,  and  execution  of  the 
blocks  in  this  list  for  each  half  time  step  is  referred  to  as  ‘execution  of  the  sort  list’.  User  output 
(UO)  blocks  arc  excluded  from  the  sort  list  because  their  output  is  calculated  when  the  associated 
user-defined  block  is  executed.  However,  when  storing  the  configuration  statements  (see  Section 


4. 1.5).  statements  with  a UO  block  type  are  stored  immediately  following  the  associated  statement 
with  a user-defined  block  type. 

The  algorithm  for  determining  the  sort  list  is  described  as  follows: 

(1)  Enter  into  the  sort  list  the  constant  (K)  blocks  in  ascending  order  of  block  number 
followed  by  the  unit  delay  (U)  blocks,  also  in  ascending  order  of  block  number. 

(2)  Select  the  lowest  numbered  unsortcd  block  and  then  go  to  (4). 

(3)  Select  the  next  lowest  unsorted  block  number. 

(4)  If  the  block  numbers  of  all  the  inputs  (i.e.  Bl,  B2,  and  B3)  have  (a)  already  been  entered 
into  the  sort  list  or  (b)  are  themselves  defined  as  integrator  (I)  or  first  order  lag  (Tl) 
blocks, t enter  this  block  into  the  sort  list  and  then  go  to  (2):  otherwise  go  to  (3). 

3.3  Model  Execution  Process 

Initially,  the  output  value  of  each  block  is  set  equal  to  the  value  of  its  first  parameter,  i.e- 
X = PI.  This  has  the  effect  of  setting  initial  conditions  for  integrator  (I),  first  order  lag  (Tl)- 
unit  delay  (U),  vacuous  (V).  and  zero  order  hold  (Z)  blocks,  and  setting  constant  (K)  blocks  to 
their  constant  value.  The  output  values  of  other  blocks  arbitrarily  set  in  this  way  are  overwritten 
on  initially  executing  the  sort  list.  The  only  difference  from  subsequent  sort  list  executions  is 
that  the  output  values  of  the  U,  V,  and  Z blocks  are  not  set  since  their  output  values  have  already 
been  set  equal  to  PI.  At  each  subsequent  half  time  step,  the  output  values  of  I and  Tl  blocks  are 
first  calculated  (sec  Section  6.2)  and  the  sort  list  is  then  executed.  Though  calculated  before 
execution  of  the  sort  list,  the  output  values  of  I and  Tl  blocks  are  not  set  until  immediately  after 
the  output  values  of  the  U blocks  are  set.  Because  the  output  value  of  the  K block  is  set  initially 
and  does  not  change,  K blocks  are  skipped  when  the  sort  list  is  executed. 

4.  BASIC  OPERATIONS  OF  SIMILATION  LANGUAGE 

This  section  explains  only  the  basic  operations  necessary  for  a user  to  be  able  to  initially 
use  the  simulation  language.  Further  operations,  which  allow  a much  greater  degree  of 
flexibility,  are  described  in  the  next  section.  The  modelling  program  BOMMP  is  first  used  to 
perform  the  simulation,  and  the  output  program  TRANS  is  then  used  to  obtain  printed  and 
plotted  output.  Both  programs  are  controlled  by  commands,  which  are  typed  following  an 
from  the  Teletype.  The  appearance  of  an  signifies  the  command  mode.  These  commands 
may  be  shortened  to  three  characters  and  are  listed  with  brief  descriptions  in  Tables  2a  and  2h. 
They  are  described  in  much  greater  detail  in  Sections  4.1.  4.2.  and  5.  For  the  modelling  program 
only,  a command  string  may  be  formed  by  a list  of  commands  separated  by  semi-colons. 

The  operating  procedures  are  best  illustrated  by  an  example.  For  convenience,  the  non- 
linear spring  problem  of  the  CSMP-9  and  II30-CSMP  manuals  is  repeated.  The  problem  may- 
be expressed  in  terms  of  the  differential  equation 

oi-  • hy  i /(y)  = 0 or  i'  = - [6\-  • f(y)]  a 

where  /(,i  ) is  a non-linear  function  defined  by  linear  interpolation  of  a set  of  coordinate  pairs 
( f,/).  A complete  description  of  the  problem  to  be  solved  is  given  in  Figure  2.  which  also  includes 
a block  diagram  representation  and  the  necessary  configuration,  parameter,  and  function 
statements. 

4.1  Running  the  Modelling  Program 

The  running  of  the  modelling  program  may  be  conveniently  divided  into  the  following 
seven  stages: 

(I)  Appropriate  relocatable  binary  files  are  loaded  into  the  computer  core  storage,  and  a 
core  image  of  the  resulting  modelling  program  is  saved. 


t The  output  value  of  integrator  and  first  order  lag  blocks  arc  known  before  the  sort  list  is 
entered  (see  Section  6.2). 
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TABI  K 2A 


Modelling  Program  C'ommand>> 

All  commands,  except  the  ‘channel’  command,  may  be  shortened  to  three  characters 


Command 

Effect 

Section(s)  for 
Reference 

CONFIGURATION 

Reads  in  configuration  statements 

3.1 

.1,  4.1.4 

DEBUG 

Provides  a complete  block  description  during 
execution 

5.5.3 

EXIT 

Returns  control  to  the  monitor 

4.1.7 

FILE 

Gives  the  status  of  the  model  input,  block  output 
and  model  output  channels 

4.1.3 

FUNCTION 

Reads  in  function  statements 

3.1 

.3.  4 14 

GOE 

Initiates  execution  of  the  model 

4.1.7 

INTEGRATION 

Reads  in  integration  control  parameters 

4.1.6 

LOOK 

Examines  block  output  at  final  time  value 

5.5.2 

MANUAL 

During  non-interactive  On-Line  processing,  returns 
the  program  to  interactive  running  and  to  the 
command  mode 

5.4 

OUTPUT 

Reads  in  output  control  parameters 

4 16 

PARAMETER 

Reads  in  parameter  statements 

3 1 

2.  4 1.4 

PROGRAM 

Executes  another  specified  program,  e g ‘TR  ANS' 

5.5.4 

RETAIN 

Enables  model  execution  to  continue  at  a later 
stage 

S S S 

RUN 

When  setting  up  a model,  initiates  the  command 
string  TITLI  :(  (fN:PAR  FUN;INT.OU  I ; 

GOE’;  if  the  Teletype  is  not  used  as  the  model 
input  channel,  the ‘TITL  ('command  is  excluded; 
where  part  of  the  string  has  already  been  initialed 
by  individual  commands,  the  remaining  ones  are 
initiated 

5 5.1 

STORE 

Stores  model  specification  statements  in  a file 

4.1 

.5,  5.1.1 

TITLE 

Reads  in  title 

4.1.4 

‘channel’ 

Used  to  change  the  status  of  the  input  and  output 
channels 

5.1 

wiih. 


(2)  Execution  of  the  program  is  initiated,  and  the  amount  of  additional  storage  is  deter- 
mined by  specifying  the  maximum  block  number  and  number  of  special  blocks  so  that 
the  appropriate  array  dimensions  may  be  expanded  to  the  minimum  size  necessary. 
The  maximum  block  number  may  exceed  the  total  number  of  blocks  used  since  not  all 
intermediate  numbers  need  be  used  as  block  identifiers. 

(1)  The  ‘status'  of  the  input  and  output  channels  may  be  examined  if  required  (see  Sections 
4.1.3  and  5.1  for  explanation). 

(4>  The  model  is  set  up  by  specifying  in  turn:  a title,  configuration  statements,  parameter 
statements,  and  function  statements  (if  any). 
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(5)  The  model  may  be  stored  for  future  use  or  reference. 

(6)  Integration  and  output  control  parameters  are  specified. 

(7)  The  model  is  run,  and  on  completion,  control  may  be  returned  to  the  monitor,  or 
modifications  may  be  made  to  the  model  by  repeating  the  appropriate  commands  so 
that  the  model  may  be  re-run. 

Stages  (3)  and  (5)  are  optional,  but  Stages  (I),  (2),  (4),  and  (6)  must  be  completed  in  order 
before  the  model  is  run  in  Stage  (7).  Stage  (3)  may  be  included  anywhere  following  Stage(2) 
while  Stage  (5)  must  be  preceded,  though  not  necessarily  immediately,  by  Stage  (4).  Stages  (3)  to 
(7)  are  controlled  by  appropriate  commands  (see  Table  2a).  Though  not  shown  here,  commands 
may  be  repeated  at  any  stage.  This  enables  alterations  or  additions  to  the  model. 

For  the  non-linear  spring  problem,  the  operating  procedures  are  described  below  for  each 
of  the  above  stages. 

4.1. 1 Stage  I — Loading  the  Program 

Aj^propriate  relocatable  binary  files,  which  include  the  files  BOMMP  Rf.L,  DUSFR  REL. 
and  any  user-defined  files,  arc  loaded  into  the  computer  core  using  the  normal  POP-IO  loading 
procedure,  and  a core  image  of  the  resulting  modelling  program  is  saved.  The  user-defined 
files  contain  user-defined  subroutines,  as  well  as  other  functions  and  subroutines  required  by  the 
user-defined  subroutines,  and  the  file  BOMMP  REL  is  the  binary  version  of  the  basic  modelling 
program.  Because  the  basic  modelling  program  contains  •C'.AEl.'  FORERAN  statements  for 
each  of  the  user-defined  subroutines  E'SERI.  USER2.  . . LiSF  RI.S  and  the  subroutines 
INTRPl,  (NTRP2,  INTRP3  used  for  user-defined  interpolation  in  the  function  block  (see 
Appendix  C),  corresponding  versions  of  these  subroutines  arc  required  when  ,.>ading  the  program. 
‘Dummy'  versions,  which  perform  no  calculations,  arc  therefore  stored  in  a file  whose  binary 
version  is  called  DUSER.REL.  By  loading  DUSl  R RE  I last  and  using  the  ‘L'  switch  (i.e. 
the  user  types  DUSER/L),  these  ‘dummy'  versions  arc  loaded  only  if  corresponding  versions 
specified  by  the  user  have  not  already  been  Kiaded.  For  the  non-lincar  spring  problem,  which 
has  no  user-defined  subroutines,  the  files  BOMMP. REI.  and  DUSFR  REL  are  loaded,  and  a 
core  image  of  the  program  is  saved  as  follows: 

R LOADER 

*BOMMP  Basic  modelling  program 

*DL'SER/L  'Dummy' versions  of  user-defined  subroutines  loaded  last 

*$  ‘Altmode'  character 

BOMMP  I2K  CORK.  194  WORDS  FREE 
LOADER  USED  14  • 4K  CORE 

EXIT 

.SAVE  BOMMP  Produces  core  image  in  file  BOMMP.SAV ; the  name  need  not 

BOMMP  SAVED  be  BOMMP 


4.1.2  Stage  2 — Execution  Initiation  and  Dimensioning  of  Expandable  Arrays 

The  execution  of  the  program  is  initiated  and  the  expandable  array  dimensions  are  set  as 
follows 


RU  BOMMP 

MAX  BLK  NO.  = 48 

NO.  OF.  / & Tl  BLKS.  U BLKS.  F BLKS  = 2,  0,  I 

II 


A 


4.1.3  Stage  3 — Status  of  Input  and  Output  Channels 

The  modelling  program  uses  three  channels  for  input  and  output.  They  are: 

(1)  an  input  channel  for  reading  in  the  model  specification  statements,  viz,  title,  configura- 
tion statements,  parameter  statements,  and  function  statements; 

(2)  a model  output  channel  for  storing  these  statements  for  future  use; 

(3)  an  output  channel  for  storing  the  block  output  (i.e.  output  value  of  each  of  the  blocks) 
in  binary  form  on  a specified  storage  device  during  execution. 

The  ‘FILE’  command  gives  the  status  of  each  of  these  channels  and  may  be  used  at  any 
stage  of  execution.  Initially,  the  channels  are  set  up  as  follows,  but  they  may  be  changed  at  any 
stage  by  using  the  ‘channel’  command  (see  Section  5.1): 

♦FIL 

MODEL  IIP  FROM  TTY: 

BLOCK  OIP  TO  LOGhOUTPT.DAT 

MODEL  OIP  TO  LOGS.  MO  DEL.  MOD 

* 


The  number  of  characters  used  for  all  file  names  (without  extension)  must  not  exceed  five. 


4.1.4  Stage  4 — Setting  up  the  Model 

The  model  is  set  up  by  specifying  in  turn  (a)  a title,  (fi)  configuration  statements,  (c)  para- 
meter statements,  and  (cl)  function  statements  (if  any).  The  appropriate  commands  ‘TITLE', 
‘CONFIGURATION',  ‘PARAMETER’,  and  ‘FUNCTION’  control  the  input  of  these 
statements  by  the  modelling  program;  the  latter  three  are  hereafter  referred  to  by  their  shortened 
form  (i.e.  ‘CON’,  ‘PAR’,  and  ‘FUN’).  Each  of  the  three  sets  of  statements  is  terminated  by 
a blank  block  number  (i.e.  an  additional  ‘carriage-return’).  The  model  for  the  non-linear  spring 
problem  may  then  be  set  up  as  shown  below,  where  typical  errors  are  included  to  illustrate  their 
correction  procedure; 

‘TIT 

TITLE  (LIMIT  60  CHRS) 

NON-LINEAR  SPRING  WITH  DASHPOT 

‘CON 

CONFIGURATIONS: 

BLK.  TYPE.  Bl,  B2,  B3 

4,  /,  17,  6;  Y DBLE  DOT 

5,  K 

6,  K 
5 

PREY  BLK  5 DELETED 

9,1,48;  Y 
17,  W,  48.  10 
48.  I.  4;  Y DOT 
10,  F,  9 

‘PAR 


Typed  in  error 

Blank  type  causes  deletion  of  statement  defined  previously 
for  block  number  5 


Blank  block  number  terminates  configuration  statements 
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r \ R \ MI  TI.RS: 


Bi.K.  /’/,  r:.  r.i 

9,  -10.0 
17,  2 
17.  2.  1 
6,  -5 

*Fl  \ 

FL  \Cr/().\S: 

Rl.K  SO.  10 

( (H) RO  FAIR: 

-10,  - KM) 

8,  -64 
-6,  -.16 
-6.  -.16 

COORD  PAIR  C/>.()()(IOt:  ■ DO.  3.6II()(IF  rOI)  DELFTFD 

-4.  -16 

0,  0 

-2,  -4  Coordinate  pairs  may  be  in  any  order 

2.  4 
8.  64 

5,  .10  T\ped  in  error 

5 Blank  ordinate  causes  deletion  without  replacement 

COORD  PAIR  ( 5.00IKIF  ■ 00.  3.00(H)E  ■ 01)  DELETED 

4.  16 

10,  100 

6.  .16 

Blank  abscissa  terminates  coordinate  pairs 

m K SO.  Blank  block  number  terminates  function  statements 

MODEL  COMPLETE 


Second  parameter  omitted  in  error 
Repeated  block  number  causes  replacement  of  parameters 
defined  previously  for  block  number  17 
Blank  block  number  terminates  parameter  statements 


i 


Coordinate  pairs  are  read  in  using  the  format  statement 
FORMAT  (21  ) 


T\ped  in  error 

Repeated  abscissa  causes  replacement 


The  message 


MODEL  COMPLETE 


(see  above  example)  is  typed  by  the  program  only  when  coordinate  pairs  for  all  the  function 
blocks  have  been  specilied.  If  there  are  no  function  blocks,  the  message  is  typed  following  the 
completion  of  the  parameter  statements.  Even  if  there  are  no  parameters  to  be  specified,  the 
PAR'  command  must  still  be  used  following  the  'CON'  command. 


4.1.5  .Stage  5 — Storing  the  Model 

By  using  the  '.STORE'  command,  the  current  model  specification  statements  may  be  stored 
for  future  use  in  a file  which  is  referred  to  as  the  model  output  channel.  For  the  present  example, 
the  complete  model  is  stored  in  'I.OG.l:MODEL.MOD'  (refer  to  initial  status  of  the  model 
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output  channel  in  section  4.1.3),  but  each  set  of  statements  may  be  stored  in  separate  model 
output  channels;  however,  the  model  output  channel  must  be  redefined  prior  to  each  set  being 
stored  (see  Section  5.1.1).  The  title  is  included  in  the  set  of  configuration  statements.  Comments 
follot^ing  the  optional  statement  terminator  are  ignored  and  hence  are  not  stored  (see 
Section  5.3  if  comments  are  required  to  he  stored).  The  appropriate  Teletype  messages  and  res- 
ponses here  are: 

•STO 

COW  PAR.  FL'W  OR  ALL  ALL 


The  configuration,  parameter,  and  function  statements  are  each  stored  in  ascending  order 
of  block  number.  However,  by  typing  CS  in  place  of  CON.  and  AS  in  place  of  ALL.  the  con- 
figuration statements  are  stored  in  sorted  order  of  block  number. 

Parameter  statements  that  have  a zero  value  for  PI,  P2,  and  P3  are  not  stored. 

4.1.6  Stage  6 — Specification  of  Integration  and  Output  C ontrol  Parameters 

The  commands  'INTFGRATION'  and  'OUTPUT',  which  are  hereafter  referred  to  by 
their  shortened  form  (i.e.  'INT'  and  'OUT'),  allow  the  integration  and  output  control  para- 
meters to  be  set;  the  order  of  these  commands  is  optional.  For  the  integration,  which  is  per- 
formed using  the  second  order  Runge  Kutta  (i.e.  Modified  Fuler)  method,  the  lower  and  upper 
time  limits  and  the  interval  are  specified  (each  in  the  same  unit).  For  the  output  control,  the 
required  block  numbers  are  first  listed.  They  may  be  typed  on  more  than  one  line,  each  line 
consisting  of  a maximum  of  72  characters;  two  consecutive  'carnage-returns’  terminate  the  list. 
If  all  the  block  numbers  are  required,  the  user  may  type  ALL  (which  may  be  shortened  to  A)  in 
place  of  the  individual  block  numbers,  in  which  case  a single  ‘carriage-return’  terminates  the 
list:  the  block  output  values  are  stored  in  ascending  order  of  block  number.  By  typing  AS  in 
place  of  ALL  or  A.  the  block  output  values  of  all  the  blocks  arc  stored  in  sorted  order  of  block 
number.  An  output  value  is  only  stored  if  it  changes  by  more  than  a prescribed  percentage  from 
the  previously  stored  value;  this  percentage  is  specified,  followed  by  the  interval  of  time  between 
output  values  in  the  same  unit  as  the  integration  parameters.  The  appropriate  Teletype  messages 
and  responses  for  the  present  example  are: 

MNT 

I NTPC  \ P ARAMS:  LOU  PR,  LPPPR.  INTERVAL  0.  20.0.1 
•OIT 


O P RLKS 

9,  4.  48.  10  Two  consecutive  ‘carriage-returns’  terminate  the  list 

O P PARAMS:  CHANGE  REQRD,  INTERVAL  = .01,  1 


If  the  output  interval  is  less  than  the  integration  interval,  the  output  values  are  recorded  at 
every  half  time  step  of  the  integration  (see  Sections  2.1  and  6.2).  To  execute  the  model  only  at 
the  initial  time  value,  the  lower  and  upper  integration  limits  should  both  be  set  equal  to  the 
initial  time  value;  the  integration  interval  need  not  then  be  set. 

On  using  the  ‘OUTPUT’  command,  a ‘TITl.F’  command  is  first  implied  if  a title  has  not 
already  been  specified  (i.e.  the  program  requires  the  user  to  type  a title). 
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4.1.7  Stage  7 — F.xecuting  the  Model 

I he  ‘Got’  eonimand  iniliales  cxeculion  of  ihe  coniplele  model.  Upon  completion,  the 
run  C PU  time  is  given  and  control  is  returned  to  the  command  mode.  If  required,  alterations 
may  then  be  made  to  the  model  before  rc-exceution  (see  Section  5.3).  Eventually,  control  is 
returned  to  the  monitor  by  the  ‘EXIT'  command.  Use  of  the  ‘GOE’  and  ‘EXIT’  commands 
is  illustrated  as  follows: 

•CiOE 


•*  KL'WIXG  ** 


Ki  \ cn-  Ti.\u::  i.:s  si:c. 

*KXI 


i:.\n  or  rxrcLTiox 

cn  riMC:  5.74  CL.-U’SCD  TlMl.:  I.C5H.74 
EXIT 


The  model  execution  may  be  halted  before  completion  b>  typing  in  which  case  execution 
IS  completed  up  to  the  next  time  value  for  which  the  block  values  specilied  for  output  are  required. 
The  message 

Rl  S TKR\US.\TED  li)  .1  


is  typed  by  the  program  and  control  is  returned  to  the  command  mode. 

4.2  Running  the  Output  Program 

I he  output  program  TR.ANS  is  controlled  by  commands  in  exactly  the  same  way  as  the 
modelling  program.  These  commands  arc  listed  with  brief  descript  ions  in  Table  2h.  The  program 
proviiles  Teletype  and  line  printer  output  in  either  tabular  or  graphical  form  and  incremental 
plotter  output  in  the  form  of  'strip'  plots  (of  the  form  produced  by  a multi-channel 
chart  recorder)  and  ‘overlay’  plots  (a  single  graph  for  up  to  eight  output  variables). 
Ihc  block  output  file  obtained  by  running  BOMMP.  which  then  becomes  the  input  file  for 
TR.ANS,  includes  the  time  limits  and  interval,  a title,  and  labels  for  each  output  variable.  Subject 
to  some  rest  rictions.  these  may  be  redefined,  tor  graphical  output,  convenient  scales  are  calculated 
by  TRANS  such  that  each  curve  is  fully  shown;  however,  these  scales  may  also  be  redefined  by 
the  user. 

Only  the  commands  necessary  to  produce  Teletype  tabular  and  graphical  output  (i.e. 
‘PRC'OTUMN’  and  ‘PRIM.OT')  are  described  here;  a complete  description  of  all  the  features 
of  TRANS  is  given  in  Reference  '). 

4.2.1  .Setting  up  the  Program 

The  User  initiates  execution  of  the  program  by  first  running  the  core  image  version 
TRANS. SA\'  of  the  relocatable  binary  file  TRANS  REl.,  and  then  specifying  the  filename 
(without  extension)  of  the  block  output  file  (OUTPT  for  the  non-lincar  spring  problem)  produced 
by  running  BOMMP.  The  output  program  then  types  the  title,  date  and  time  of  creation  of  the 
block  output  file,  integration  interval,  run  C PU  time,  and  time  parameters,  followed  by  an 
which  signifies  the  command  mode.  An  example  of  how  to  set  up  the  program  in  the  above 
manner  is  shown  below. 


TABLE  2B 

Output  Program  Commands 

All  commands  may  be  shortened  to  three  characters 


1 


Command 

Effect 

Section(s)  for 
Reference 

EXIT 

Returns  control  to  the  monitor 

4.2,1 

GOE 

Obtains  the  tabular  or  graphical  output  for  the  commands 
‘PRCOLUMN’,  ‘PRPLOr,  PLSTRIP',  and 

•PLOVERLAY’ 

4.2.1 

LABEL 

Redefines  labelling  information 

« 

PLOVER LAY 

Specifies  ‘overlay’  plotting  on  the  incremental  plotter 

* 

PLSTRIP 

Specifies  ‘strip’  plotting  on  the  incremental  plotter 

* 

PRCOLUMN 

Specifies  tabular  output  on  the  disk  or  Teletype  in  the 
form  of  printed  columns 

4.2.2 

PRPLOT 

Specifies  graphical  output  on  the  disk  or  Teletype 

4.2.2 

RUN 

Equivalent  to  the  commands  (in  order)  ‘PRCOLUMN', 
‘PRPLOr,  ‘PLSTRIP’,  and  ’GOE’;  if  any  form  of 
output  is  not  required,  a ’carriage-return’  is  typed  in 
place  of  the  block  numbers 

* 

SCALE 

Redefines  graphical  output  scales 

* 

TIME 

Redefines  time  parameters 

♦ 

* See  Reference  9. 


RU  TRANS 

IP  FILENAME  = OL’TPT 

NON-LINEAR  SPRING  WITH  DASHPOT 

IP  FILE  RECORDED  ON  2I-APR-76  AT  16:57 

INTEGN  !NT=  I.OOOOE-OI,  RUN  CPU  TIME  = 1.25  SEC 

TIME  FROM  O.OOOOE-OI  TO  2.nOOOE+OI  IN  STEPS  OF  ! .mWE - 00 


Although  the  filename  extension  is  not  specified,  the  extension  name  ‘OAT’  is  assumed.  If. 
in  the  above  example,  a file  named  OUTPT.DAT  is  not  found  on  the  disk,  the  message 

OUTPT.DAT  NOT  ON  DSK 

UP  FILENAME^ 
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A... 


IS  lypcJ  by  the  profiram  so  lluil  ;i  new  lilenanie  may  be  speeilied. 

A parlieiilar  type,  or  types,  of  output  may  be  speeilied  by  using  the  appropriate  commands 
listed  in  Table  2b:  to  complete  the  output  process,  the  'GOI-'  command  should  follow,  after 
which  control  is  returned  to  the  command  mode.  Repetition  of  the  appropriate  commands  for 
each  type  of  output  before  the  ’GOfi’  command,  deletes  the  ciTect  of  the  previous  use  of  the 
command ; repetition  after  the ‘GOIi’ command,  results  in  additional  output.  Control  is  returned 
to  the  monitor  by  the  '[{XIT'  command. 

^ 4.2.2  I'eletype  Output 

The  commands  'PRCOITIMN'  and  ‘PRf’l.Or  are  used  to  obtain  tabular  and  graphical 
output  respectively  on  either  the  line  printer  or  Teletype.  Tor  each  command  the  required  block 
numbers  are  first  listed.  They  may  be  typed  on  mcire  than  one  line,  each  line  consisting  of  a 
masimum  of  11  characters,  and  the  list  is  terminated  by  two  consecutive  ‘carriage-returns'.  If 
all  the  bloek  numbers  are  required,  the  user  may  type  A in  place  of  the  individual  block  numbers, 
in  which  case  a single  'carriage-return'  terminates  the  list.  When  specifying  blocks  by  their 
individual  block  numbers,  a maximum  of  a hundred  blocks  may  be  processed  each  time  the 
'CiOl  ' command  is  used:  however,  all  blocks  are  processed  when  the  user  types  instead,  even 
if  there  are  more  than  a hundred  blocks,  following  the  Telety  pe  message 

IS  O />  TO  TTY  RTQRP: 

typed  by  the  program,  the  user  types  \ if  the  output  is  tii  be  printed  directly  on  the  Teletype,  or 
any  other  character  (including  a 'carriage-return')  if  the  output  is  to  be  stored  in  a disk  file  so 
that  the  resulting  tile  can  subsequently  be  printed  on  the  line  printer  (see  Reference  9),  This  latter 
form  of  output  contains  more  information  than  the  Teletype  output;  blocks  are  identified  by 
their  number  rather  than  by  a label  on  the  Teletype  output.  On  completing  the  Teletype  responses 
to  either  the  'PRtOI  LiMN'  or  'I’RPI  Of  commands,  the  program  returns  control  to  the 
command  mode.  The  'GOT.'  command  is  then  nsevl  to  obtain  the  Teletype  output.  It  should 
be  noted  that  both  tabular  and  graphical  output  cannot  be  obtained  on  the  Teletype  by  using  a 
single 'GOT' command  following  the  'PRC'OflJMN' and  'PRPI  OT' commands.  Convenient 
plotting  scales  are  calculated  by  the  output  program  ;ind  are  shown  at  the  head  of  each  Teletype 
plot . 

Ihe  following  example  shows  how  both  tabular  and  graphical  Teletype  output  may  be 
I'htained  for  the  non-linear  spring  problem: 

*I>R( 

TR/MISa  l\  COl.LMXS: 

HISS 

A fist  terminated  by  a single  'carriage-return' 

IS  O T TO  TTY  REQRD:  \ 

*(;oK 

**  R(  \M\(i  ** 

(See  Sec(ion  I of  Appendix  I)  for  Teletype  output) 

♦PRP 

I’RIMTR  ri.OTS: 


fist  terminated  by  two  consecutive  ‘carriage-returns’ 


RTKS 

9.  4 
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>( 


IS  O P TO  TTY  REQRD:  Y 
*GOE 

**  RUNMING  ** 

(See  Section  2 of  Appendix  D for  Teletype  output) 

*EXI 

END  OE  EXECUTION 

CPU  TIME:  2.27  ELAPSED  TIME:  9:9.H2 

EXIT 


5.  FURTHER  OPERATIONS  OF  MODELLING  PROGRAM 

In  the  previous  section,  the  basic  operations  necessary  to  run  both  the  modelling  program 
and  output  program  were  described.  In  this  section,  further  operations  of  the  modelling  program 
are  described  that  allow  much  greater  flexibility.  In  particular,  facilities  are  described  that  allow 
repeated  alteration  to  the  model  and  repeated  execution  of  the  model  without  exiting  from  the 
modelling  program. 

5.1  Changing  the  Status  of  Input  and  Output  Channels 

The  three  channels  used  for  input  and  output  are  described  in  Section  4.1.3,  where  the  initial 
status  of  each  channel  is  given  by  the  'FILE'  command.  If  the  user  wishes  to  re-run  a model, 
perhaps  after  making  some  alterations,  it  may  be  desirable  to  alter  the  status  of  one  or  more  of 
these  channels.  The  logical  unit  and  filename  for  a particular  channel  is  altered  by  the  program 
command  (termed  a ‘channel’  command) 

I.OGn:  “filename"-  “channel" 

where 

(1)  n is  an  integer.  I to  12,  representing  the  logical  unit  number;  the  following  list  gives 
those  that  are  reserved,  LOGS  being  permanently  reserved: 

l.OGI  - block  output  channel  (initially); 

LOGS  -model  output  channel  (initially) 

LOGS  — output  file  named  BOMMP.OUT  which  contains  the  Teletype  responses 
typed  by  the  user  for  subsequently  running  the  modelling  program  non- 
interactively  (see  Section  S.4); 

LOG  1 2— output  file  named  DEBUG. OUT  which  contains  a complete  block 
description  (i.e.  values  for  each  block  of  X,  XI,  X2,  X3,  M,  PI,  P2,  and 
P3)  for  each  half  time  step  during  execution;  the  file  is  only  opened,  and 
logical  unit  only  reserved,  following  an  initial  ‘DEBUG'  command 
(sec  Section  S.S.3)  when  the  output  is  not  required  on  the  Teletype. 

(2)  “filename"  is  a filename  without  an  extension  and  with  a maximum  of  five  characters. 

(3)  “channel”,  which  is  one  of  the  following  list,  represents  a particular  input  or  output 
channel : 

I'P—  input  of  model  specification  statements  (to  run  the  model); 

MODEL— output  of  model  specification  statements  (to  store  the  model); 

O/P-  block  output. 

The  above  channel  descriptors  may  be  shortened  to  one  character  (e.g.  ‘I/P’ may  be  short- 
ened to  ‘I’). 

The  Teletype  may  be  used  as  the  model  input  or  model  output  channel,  in  which  case  the 
‘channel’  command  takes  the  form 


The  Teletype  is  initially  used  as  the  model  input  channel.  The  model  specification  statements 
entered  in  this  way  are  separate  from  the  commands  and  control  parameters  read  in  using  the 
Teletype.  The  status  of  the  model  output  and  block  output  channels  remains  the  same  unless 
changed  by  a further  channel  command.  However,  if  the  status  of  the  model  input  channel  has 
been  changed  from  the  Teletype  to  reading  from  a file,  it  reverts  to  the  Teletype  once  input  is 
completed  (see  Examples  I and  2 in  Section  5.1.2).  If  further  specification  statements  are  required 
to  be  read  in  from  another  file,  the  status  of  the  model  input  channel  would  have  to  be  changed 
from  the  Teletype  to  the  new  tile. 

Whenever  the  status  of  a channel  is  changed,  either  by  a ‘channel’  command  or  upon  rever- 
sion to  the  Teletype  for  the  model  input  channel,  the  status  of  the  appropriate  channel  is  typed  on 
the  Teletype  by  the  program  as  shown  in  the  examples  in  Sections  5.1.1  to  5.1.3. 

The  use  of  dilfcrent  logical  units  enables  tiles  to  be  read  from,  or  written  onto,  a specific 
device  other  than  the  disk  (i.e.  DF.Ctape  or  magnetic  tape).  To  use  this  facility,  the  'ASSIGN' 
monitor  command  is  used;  for  example,  to  assign  DECtape  OTA3  to  LOG  I.  the  user  types 
AS  DTA3  1 to  the  monitor  prior  to  running  the  modelling  program. 

5.1.1  Model  Output 

In  the  non-linear  spring  problem  of  Section  4,  by  using  the  'STORE'  command,  the  model 
specification  statements  were  stored  in  ‘LOG3:  MODEL. MOD'  (see  Sections  4.1,3  and  4.1.5). 
The  'STORE'  command  allows  any  one.  as  well  as  all,  of  the  separate  groups  of  specification 
statements  (i.e.  configuration,  parameter,  and  function  statements)  to  be  stored  in  the  file  defined 
as  the  model  output  channel.  The  three  examples  below  show  how  (I)  all  the  specification  state- 
ments may  be  stored  in  'LOG6:SPRlNG.MOD',  (2)  the  configuration,  parameter,  and  function 
statements  may  be  stored  separately  in  'LOG7:SPCON.MOD',  'LOG8:SPDAT.MOD', 
and  ‘LOG7:SPFUN.MOD‘  respectively,  and  (3)  all  the  specification  statements  may  be  listed 
on  the  Teletype  by  using  the  Teletype  as  the  model  output  channel  (headings  are  deleted). 

Examph’  I 

*1.0G6:SPRNG-  M 

MO  DEI.  OP  TO  I.OG6:SPRSG.MOn 

*.STO 

COS.  PAR.  TVS.  OR  ALL  ALL 


Example  2 

•LOr.7:SPCON-  M 

MODEL  O P TO  LOG7:SPCON.MOD 

*.STO 

COS.  PAR.  PCS.  OR  ALL  .CON 
‘LOG8.SPDAT.  .VI 

MODEL  O P TO  l.OGS.SPDAT.MOD 

*STO 

COS,  PAR.  EUS.  OR  ALL  PAR 
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LOG7;SPFl  N-  M 


1 


MODEL  O P TO  LOOT: SPED. MOD 


C'O.V,  PAR.  EU.X.  OR  ALL  Fl'N 


Example  3 

*TTV:-  M 

MODEL  O P TO  TTY: 

•sro 

cos.  P.iR.  El  S.  OR  ALL  . AS  Configuration  statements  in  sorted  order  (see 

Section  4. 1 .5) 

.\0.\-LI.\EAR  .SPRLSC  W ITH  DASHPOT 
6 K 


d 

/ 4S 

y 

1(1 

E 9 

17 

IT  43 

10 

4 

1 /7 

6 

;Y  DBLE  DOT 

4H 

/ 4 

. }'  DOT 

6 

5.0000E  ■ 

IK) 

9 

I .OOOOE  ■ 

01 

17 

2.0000E  ■ 

00 

I.OOOOE  ■ 00 

10 

I.OOOOE 

■ 01 

I.OOIKIE 

02 

H.moOE 

\ 00 

6.4000E 

01 

6. 0000 E 

. 00 

J.6000E 

01 

4.0000E 

■ 00 

I.6000E 

01 

- 2. 0000 E 

■ 00 

4.0000E  • 

00 

O.OOOOE 

01 

O.OOOOE 

01 

2.0000E 

■ 00 

4.0000E 

00 

4.0000E 

■ 00 

L6000E  • 

01 

O.OOOOE 

■ 00 

3.6000E 

01 

3.0000E 

■ 00 

6.4000 E • 

01 

I.OOOOE 

■ 01 

I.OOOOE  ■ 

02 

* 


Usually,  when  reading  from  a file  used  as  the  model  input  channel,  all  the  statements  would 
be  read  in.  However,  it  is  possible  that  the  different  types  of  statements  may  need  to  be  read  in 
Irom  dilTerent  sources  (i.e.  Teletype  and  various  files)  without  reading  in  all  the  information  stored 
in  a particular  file.  For  example,  two  files  may  each  contain  a set  of  configuration,  parameter, 
and  function  statements.  The  user  may  wish  to  read  the  configuration  statements  from  one  file 
(File  A say),  fhe  parameter  statements  from  the  Teletype,  and  the  function  statements  from  the 
other  tile  (File  B say).  The  procedure  would  be  as  follows: 

(1)  Ifse  tfic  ‘channel’  command  so  that  File  A is  used  as  the  model  input  channel. 

(2)  Type  CON  to  read  in  the  configuration  statements. 
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(3)  Use  the  'channer  command  so  that  the  Teletype  is  used  as  the  model  input  channel; 
the  parameter  and  function  statements  are  therefore  not  read  in  from  File  A. 

(4)  Type  PAR  followed  by  the  parameter  statements. 

(5)  Use  the  ‘channel'  command  so  that  File  B is  used  as  the  model  input  channel, 

(6)  Type  FUN  to  read  in  the  function  statements;  the  configuration  and  parameter  state- 
ments are  therefore  skipped  in  File  B. 

Although  the  user  is  free  to  redefine  the  data  source  for  each  type  of  statement,  it  should  be 
remembered  that  when  setting  up  the  modrl,  the  commands  ‘CON’,  ‘PAR’,  and  ‘FUN’  (if 
any  functions)  must  be  typed  in  order;  the  ‘TITLE'  command  would  also  be  required  if  reading  the 
configuration  statements  from  the  Teletype.  Until  the  model  is  complete,  modifications  should  not 
generally  be  attempted  by  repeating  the  commands  (e.xcept  for  the  ‘TITLE’  command).  The  model 
input  channel  only  reverts  to  the  Teletype  on  reacning  •ihee"d  of  a file.  If  a file  is  not  read  to  its  end 
therefore,  as  in  the  case  of  File  A in  the  above  e.xample,  a ‘chain, t,'  command  must  be  used 
to  either  revert  to  the  Teletype  or  change  to  another  tile.  On  mai  ng  an  error  in  the  above 
procedures,  an  appropriate  diagnostic  error  message  (Appendix  B)  is  given,  and  control  is  returned 
to  the  command  mode  with  the  model  unaffected. 

5.1.2  .Model  Input 

Having  already  stored  the  specification  statements  during  a previous  run  of  the  modelling 
program,  the  user  may  wish  to  set  up  the  model  using  these  stored  statements.  Alterations  may 
then  be  made  after  all  the  statements  have  been  read  in  (see  Section  5.3)  from  the  current  input 
file.  Because  the  title  is  included  in  the  configuration  statements  when  using  the  ‘STORE’ 
command,  the  ‘TITLE’  command  should  not  be  used  unless  the  titie  is  to  be  modified  (see 
Section  5.3),  in  which  case  it  should  be  used  after  the  configuration  statements  have  been  read  in. 
If  it  is  used  before,  the  title  included  in  the  configuration  statements  will  overwrite  the  intended 
title.  Consider  Examples  I and  2 given  in  Section  5.1.1,  where  the  same  specification  statements 
are  stored  in  two  dilferent  ways.  The  method  of  reading  these  statements  for  both  examples  is: 

Example  I 

*LOG2:SPRNG-  I 

MODEL  IIP  FROM  l.OG2:SPR\G.MOD 

*CON 

CONFIGVRATIOSS  : 


SON-LINEAR  SPRING  WITH  DASH  POT 


BLK 

TYPE  BI  B2 

BJ 

4 

IT  6 

, T DBLE  DOT 

6 

K 

9 

1 4H 

y 

10 

F_  9 

17 

w'  48  10 

4,S 

! 4 

, >•  DOT 

*PAR 

PARAMETERS  : 

BLK 

P! 

P2  P3 

6 

-5. 00(H) E • 00 

9 

-UHHH)E  ■ 01 

17 

2.(HHH)E  f (H)  I.OOOOE  ■ 00 

Kl  N 


n SCTIOSS  : 

Bl  .K  .\  0.  10 

COORD  PAIRS  : 


-LOOOOE  ■ 01 

N. OOOOE  ■ 00 
6. 0000 1:  00 
4. 0000 F.  ■ 00 

:.ooooF  ■ 00 

O. OOOOF  01 
2.0000F  ■ 00 
4.0000F  ■ 00 
O.OOOOF  00 
S.OOOOF  00 
I .OOOOF  ■ 01 


I OOOOF  ■ 02 
6. 4000  F 01 
3.6000F  ■ 01 
I. OOOOF  01 
4.0000F  ■ 00 
O.OOOOF  01 
4.0000F  00 
I.OOOOF  01 
S.OOOOF  ■ 01 
0.4000F  ■ 01 
I.OOOOF  ■ 02 


MODFL  IP  FROM  L0(J2:SPRSG.M0D  COMPFFTFD 

MODFL  IP  FROM  TTY:  Status  of  model  input  channel  reverts  to  Teletype 


MODFL  COMPLFTF 


Example  2 

»LOG6:SPCON  I 

MODEL  IP  FROM  LOOO  SPCON  MO D 

*co\ 

CO Miau RATIO \S  : 

i\0.\'-LI.\FAR  SPRI.Xa  WITH  DASIIPOT 


I.F 

TYPE 

B! 

B2  BS 

4 

0 

K 

17 

6 :Y  DBLF  DOT 

9 

1 

4S 

) 

10 

1 

9 

17 

H 

4S 

10 

4,S 

1 

4 

: Y DOT 

MODFL  IP  FROM  LOOO.SPCOS'.MO D COMPLETED 
MODFL  / P FROM  TTY: 

»I,(K;6:SPI)AT-  I 

MODFL  IP  FROM  LO(,0:SPD.4T.MO D 

*PAR 
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PARAMETERS  : 


BEK  PI  P2  P3 

6 S.OOOOE  t 00 
9 -I.OOOOE  0! 

17  2. 0000 E ■ 00  I.OOOOE  00 

MODEL  IP  ERO.U  LOG6:SP D.iT.MO D COMPLETED 

MODEL  ! P EROM  TTY: 

M,OC;6.SPHN-  I 

MODEL  IP  EROM  L006:SPEU\..MO D 

*FLN 

I i WCTlOyS  : 

BI.K  \0.  10 


COORD  PAIRS  : 


I.OOOOE  ■ 01 
R.OOOOE  ■ 00 
- 6. 0000 E ■ 00 
-4.0000E  ■ 00 
-2. 0000 E ■ 00 
O.OOOOE-OI 
2.0000E  ■ 00 
4.0000E  ■ 00 
6. 0000 E ‘ 00 

H. OOOOE  i 00 

I. OOOOE  - 01 


-I.OOOOE  ■ 02 
-6. 4000 E ■ 01 

- 3.6000E  • 01 
-I.6000E  ■ 01 

- 4.0000E  ; 00 
O.OOOOE-OI 
4. 0000 E > 00 
I.6000EA  01 
3.6000E  ■ 01 
6.4000E  ■ 01 
I.OOOOE  02 


MODEL  IP  EROM  LOG6:SPEC N MO D COMPLETED 


MODEL  IP  EROM  TTY: 
MODEL  COMPLETE 


5.1..^  Block  Output 

While  running  the  modelling  program,  the  user  may  wish  to  make  a number  of  separate 
alterations  either  to  the  specitication  statements  or  control  parameters,  and  to  re-execute  the  up- 
dated model  cadi  time.  To  preserve  the  block  output  for  each  run,  the  block  output  channel  may 
be  altered  as  in  the  following  example: 

*I  (H;I:01TI-  o 

BLOC  K (>iP  TO  LOG!  OCT!  DAT 
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i 

1 


5.2  Kile  Safety 

5.2.1  Output  Files 

ir  the  user  attempts  to  overwrite  a model  or  block  output  tile  that  is  already  in  existence  on 
the  disk,  then  the  I'ollowing  message  is  typed  on  the  Teletype  by  the  program: 

BLOCK  MODELf  O P TO  TOCn, "filename". "extension" 

MILL  OVERli'RITE  “filename". "extension"  ON  DSK 
OK 

If  the  file  is  to  be  over^vntten  or  written  onto  a device  other  than  the  disk,  having  previously 
assigned  the  device  to  the  appropriate  logical  unit,  the  user  types  V and  the  program  will  continue. 
Otherwise,  the  user  types  any  other  character,  including  a ‘carriage-return'.  The  message 

THEN  CHANGE  EILENAME 

is  then  typed  by  the  program  and  control  is  returned  to  the  command  mode.  The  status  of  the 
model  or  block  output  channel  should  then  be  altered  by  using  the  appropriate  ‘channel’ 
command  (see  Section  5.1)  before  continuing.  Block  output  and  model  output  files  are  opened 
automatically  following  the  ‘GOE‘  and  ‘STORE"  commands  respectively. 

Following  a system  failure  during  Batch  processing,  it  may  be  possible  that  output  files 
opened  during  the  ‘failed  run'  would  need  to  be  overwritten  when  the  run  is  repeated,  thus 
requiring  interaction  by  the  operator.  To  avoid  this,  when  running  the  modelling  program  non- 
interactively  (see  Section  5.4),  an  output  file  will  be  overwritten  automatically;  in  place  of  the 
above  messages,  the  following  message  is  typed  by  the  program: 

HAVE  OVERH'RITTEN  FILE  "filename". "extension"  ON  DSK 

5.2.2  Input  File 

Immediately  following  a channel  command  that  changes  the  status  of  the  model  input 
channel,  the  modelling  program  searches  for  the  appropriate  filename  on  the  disk.  If  it  cannot 
be  found,  the  following  message  is  typed  by  the  program: 

MODEL  EP  FROM  LOGnCMename" .MO D 

“filename".. t/OD  NOT  ON  DSK 

IS  IT  ELSEM  HERE  : 

If  the  file  is  to  be  read  from  a device  other  than  the  disk  (see  Section  5.1),  the  user  types  Y and 
the  program  will  continue.  Corresponding  to  the  procedure  with  output  files,  by  typing  any 
other  character,  the  message 

THEN  CHANGE  FILENAME 

is  typed  by  the  program  and  control  returns  to  the  command  mode  so  that  the  input  file  status 
may  be  altered  before  continuing. 

5.it  Modifying  the  Model 

The  specification  statements  stored  in  a file  may  be  altered  or  added  to  using  a text  editing 
program  such  as  ‘TF.CO'  or  ‘EDIT’  before  running  the  modelling  program.  If  it  is  desired 


' Whichever  is  appropriate. 
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to  retain  the  comments  following  the  dollar  symbol  or  to  retain  an  unsorted  order  of  statements. 
It  is  necessary  to  make  alterations  in  this  way.  Alternatively,  specification  statements  may  be 
altered  or  added  to  by  repeating  the  appropriate  command  (i.e.  ‘TITLE’,  ‘CON’,  ‘PAR’, 
or  'EUN')  at  any  stage  after  the  model  has  been  set  up;  the  modifications  would  usually  be  read 
in  from  the  Teletype,  but  they  may  be  read  in  from  a file  that  is  being  used  as  the  model  input 
channel  (see  Sections  5. 1 and  5. 1 .2).  Alterations  may  also  be  made  to  any  of  the  control  parameters 
by  repeating  the  appropriate  command  (i.e.  ‘INT’  or  ‘OUT’),  or  to  the  status  of  the  input  and 
output  channels.  Following  these  alterations,  the  model  is  re-executed  using  the  ‘GOE’ 
command. 

5.4  Non-Interactivc  Running  and  Batch  Processing 

To  run  the  modelling  program  in  a non-interactive  way,  all  the  Teletype  input  that  is  typed 
by  the  user  when  running  interactively  must  be  stored  in  a file  named  BOMMP.IN.  The 
modelling  program  first  searches  for  a file  with  this  name  on  the  disk;  if  it  is  found,  then  the 
usual  Teletype  input  is  read  from  this  file  and  most  messages  usually  typed  by  the  program  are 
suppressed.  This  facility  has  been  introduced  primarily  for  Batch  processing,  but  may  also  be 
used  during  On-Line  processing  to  avoid  tedious  re-typing  of  the  Teletype  input  by  the  user 
when  It  is  desired  to  run  a model  that,  except  for  minor  modifications,  has  been  run  previously. 
When  using  this  facility  during  On-Line  processing,  the  messages  typed  by  the  program  are 
shown  in  the  following  example: 

Rl  BOMMP 

\OS-UNEAR  SPRING  WITH  DASH  POT 
**  RUNNING  ** 

RUN  CPU  TIME  : 1.21  SEC. 


END  OE  EXECUTION 

CPU  TIME:  2.62  ELAPSED  TIME:  3.10 

EXIT 


The  most  convenient  way  to  create  a file  containing  the  Teletype  responses  is  to  perform  a 
'dummy'  run  using  the  modelling  program  interactively  with  all  quantities  set  to  their  desired 
value  except  for  the  upper  time  limit  of  integration,  which  is  set  much  smaller  than  desired  so 
that  a minimum  of  CPU  time  is  used.  The  Teletype  responses  typed  by  the  user  are  stored  by  the  , 

program  in  ’LOG5:BOMMP.OUT‘.  The  file  BOMMP.OUT  may  then  be  renamed 
BOMMP.IN  and  the  upper  timelimit  of  integration  changed  to  the  desired  value  using  a text  \ 

editing  program.  This  method  of  creating  the  file  provides  the  user  with  a check  on  the  file  format 
so  that  he  can  be  sure  the  modelling  program  will  run  properly  when  used  non-interactively. 

To  obtain  the  advantages  of  both  non-interactive  and  interactive  running  during  On-Line 
processing,  the  ‘MANUAL'  command  has  been  introduced,  which  returns  the  program  to 
interactive  running  and  to  the  command  mode.  This  enables  alterations  to  be  made  interactively  j 

to  a model  that  has  been  read  in  non-interactively  and  allows  repeated  model  executions  follow- 
ing further  model  alterations.  Even  if  no  alterations  are  envisaged,  replacement  of  both  the 
‘GOE’  and  ‘EXIT'  commands  by  a single  'MANUAL'  command  is  advised  and  will  require  j 

only  the  'GOE'  and  'EXIT'  commands  to  be  typed  interactively  by  the  user.  If  the  ‘MANUAL' 
c’ommand  is  typed  while  running  interactively,  e.g.  when  creating  a file  to  be  renamed 
BOMMP.IN  by  performing  a dummy  run  (see  above),  control  is  returned  to  the  command 
mode  (ic.  in  elTcct.  nothing  happens  except  that  the  command  is  written  into  the  file  j 

BOMMP  OUT;  see  Section  .5.1).  On  returning  the  program  to  interactive  running  by  using  the  | 


2.5 


'MANUAL'  commatKl,  further  Teletype  input  responses  typed  by  the  user  are  stored  in  the 
tile  LROMAN.OLJT  (on  LOG5),  and  not  in  BOMMP.OUT 


I 


5.5  Miscellaneous  Commands 

5.5.1  Rl  N 

When  setting  up  a nit>del.  the  'RUN'  command  may  be  used  to  initiate  the  command 
string  ‘TlTl-L  ;CON  :PAR  ;l  LIN  :l NT: OLiTiCiOl  ' If  ihe  Teletype  is  not  used  as  the 
model  input  channel,  the  TITI  L'  command  is  escludcd  because  the  title  is  included  with  the 
configuration  statements  when  previously  stored  using  the  'STORE  ' command.  Where  part 
of  the  string  has  already  been  initialed  by  individual  commands,  the  remaining  ones  are  initiated. 
If  required,  the  'SIORI  ' command  may  then  be  used  following  the  eveculion  Having  set  up  a 
model  and  executed  it  once,  the  'RUN'  command  is  equivalent  to  the  'OOE:'  command  and 
could  be  used  following  aliertiiions  to  the  model  by  other  commands. 

If  a file  is  used  as  the  model  input  channel,  the  configuration  (including  the  title),  parameter. 
;ind  function  statements  are  all  typed  by  Ihe  program  in  response  to  the  'RUN'  command.  The 
felety  pe  output  of  these  statements  may  be  terminated  by  typing  "Control  O".  If  all  the  model 
specilieation  statements  are  not  stored  in  the  file  used  as  the  model  input  channel,  a message 
indicating  that  a particular  type  of  specification  statement  cannot  be  found  is  typed  by  the 
program  and  control  is  returned  to  the  command  mode.  1 he  model  input  channel  should  then 
be  altered  the  appropriate  number  of  limes  so  that  the  remaining  specification  statements  may 
be  read  in  on  repeating  the  'RL'N'  camimand. 

5.5.2  I.OOK 

l ollowing  the  execution  of  a model,  the  output  value  of  any  block  may  be  examined  by 
using  the  'LOOK'  command.  Tor  the  non-linear  spring  problem,  an  example  of  its  use  is  as 
follows,  control  being  returned  to  Ihe  command  mode  on  typing  a blank,  or  zero,  block  number: 

*l;ok 

*•  RC.VA7.VG  ** 

R'  V cm  n\ti:  . i.:j  six 

*ioo 

lil.K  4 

()  r i.oisii:  III 

R/R  <) 

O R 1. 621101:  01 

HIK  10 

o /’  j.2.mi:  01 


5.5.^  DKBl  (. 


I 


>< 


lor  eadi  block,  ihc  'D!  hL^Ci'  command  pro\idc>  the  block  output  value  X.  the  values 
\ I . \2,  and  X V the  mteeer  parameter  M.  and  the  parameter  values  F’l,  P2,  and  F’3  (for  notation, 
see  Section  .VI).  I he  above  values  are  printed  for  each  halftime  value  during  e.xecution  of  the 
sort  list  In  response  to  prompts  on  the  Teletype,  the  user  specifies  (u)  the  time  interval  during 
which  the  debugging  iiutpiit  is  required  and  (M  whether  the  output  is  required  directly  on  the 
Feletvpe  If  the  resp»mse  to  (/>)  is  anything  other  than  ^ . then  the  file  named  DEBUG. OUT 
IS  opened  on  I ()Cil2  if  it  has  not  already  been  opened  by  a previous  ‘FDEBUG’  command.  As 
with  the  input  ,ind  output  channels,  the  file  DF  HUG. OUT  may  be  written  onto  a specific 
device  other  th.in  the  disk  by  assigning  the  appri'priate  device  to  F.OGI2  (see  Section  5.1). 

\fter  the  model  has  been  evecuted.  debugging  for  further  runs  is  cancelled.  The  'DEBUG' 
ci'mm.ind  must  therefore  be  repeated  prior  to  ,iny  further  runs  if  debugging  output  is  again 
required  fo  obtain  debugging  output  at  one  time  value  only,  the  lower  and  upper  limits  should 
be  made  the  s.ime  I or  the  iw'n-lmear  spring  problem,  the  use  of  the  'DEBUG'  command  is 
illustr,ited  by  the  following  evample: 

»l)f  B 

//'//  liuirs  LOU  I K.  l /’/’I. R 0.2.0..^ 

IS  O r lO  in  Rl.QRH  \ 

*(,()K 

Rl  \ \/\(, 

(See  Appendix  F for  Teletype  output) 

Ri  A cm  rnii:  . 2.27  si:c 

* 


5.5.4  F*R()(.K  \M 

1 he  'F’ROGRANT  cimim.ind  enables  another  program  (e  g.  FRANS)  to  be  run  without 
specifically  terminating  the  execution  of  BOMMFV  Once  execution  of  the  new  program  is  com- 
pleted. control  IS  returned  to  the  monitor,  not  BOMMFV  Following  the  'PROGRAM' 
ci'nimand.  the  message 

RROCRW!  S I Mi- 
ls typed  on  the  Teletype  by  the  program.  The  user  responds  by  typing  the  name  (maximum  of 
live  characters)  ot  atu'ther  program.  F ollowing  a 'carriage-return',  the  new  program  is  executed. 

5.5.5  KETAFN 

I or  e.ich  block  with  type  I.  T I.  U.  and  /.  (i  e.  those  requiring  specified  initial  conditions), 
the  'RF  F AIN'  command  resets  the  initial  condition  ti  e the  first  parameter)  to  the  last  calcu- 
lated output  value.  When  used  fcdlowing  execution  of  the  model,  the  'RETAIN'  command 
enables  execution  to  be  coniinueil  ;it  a later  stage  during  either  the  same  program  run  or  a sub- 
sequent one.  In  the  latter  case,  the  'STORF  ' command  would  have  to  be  used  to  retain  the 
modilied  initi.il  cimditions  in  a file  for  future  use. 

6.  lOOl'  BREAKINt.  Bl  (K  kS 

When  solving  dilferential  or  dilFerence  equations,  it  is  necessary  to  break  the  sort  loop 
that  occurs.  Blocks  whose  output  values  are  known  at  the  Fvcginning  of  each  half  time  step 
enable  this  to  be  achieved  F I'r  diflercntial  equations,  integrator  and  First  order  lag  blocks  are 
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used,  vshilc  for  dilVerence  equations,  unit  delay  blocks  are  used.  Only  brief  descriptions  of  these 
blocks  are  gisen  in  Table  I . further  details  are  given  below. 

Ihe  vacuous  (V)  block  is  used  to  break  the  sort  loop  that  occurs  when  solving  an  implicit 
equation.  However,  this  block  can  only  be  used  in  conjunction  with  the  wye  (Y)  block;  its  use 
IS  therefore  described  as  a special  facility  in  Section  8.1. 

6.1  I'nit  Delay- 

Unit  delay  (U)  blocks  are  entered  in  order  of  ascending  block  number  at  the  beginning  of 
the  sort  list  following  the  constant  (K)  blocks  and  arc  the  first  to  be  executed  at  each  step  (see 
Section  .1.2).  This  is  because  the  input  signal  to  a U block  is  calculated  at  the  previous  half  time 
step  (see  definition  in  Table  I)  so  that  it  is  not  necessary  for  the  input  Bl  to  have  already  been 
sorted.  Hence,  provided  the  input  of  the  U block  is  not  itself  another  U block,  the  output  value 
X of  the  U block  is  delined  simply  as  the  output  value  XI  of  the  input  block  Bl.  The  procedure  is 
more  complex  when  the  input  to  a U block  is  itself  a U block;  this  occurs  when  the  output  value 
of  a block  at  the  prev  ious  complete  time  step  is  required.  The  validity  of  the  above  definition  of 
the  function  performed  by  a U block  then  depends  on  the  order  of  entry  of  the  U blocks  in  the 
sort  list.  To  avoid  this  dilTiculty.  the  U blocks  are  therefore  executed  in  the  following  two  stages, 
with  the  first  stage  being  completed  for  all  the  U blocks  before  commencing  the  second  stage 
immediately  after: 

(a)  Sla^f  I 

t When  I f,,'  ^ I’' 

I When  I > t„, 

X XI  if  B I is  not  a U block 

X P2  if  Bl  is  a U block 

(Both  XI  and  P2  are  last  calculated  at  previous  half  time  step) 

(/))  .S'rugc  2 (for  all  values  of  i) 

No  further  calculation  if  Bl  is  not  a U block 

P2  XI  if  B I is  a U block  (P2  is  then  used  in  Stage  I at  next  half  time  step) 

Consider  a function  r.  which  is  calculated  in  a user-defined  block,  U.1  say,  and  assume  a 
unit  time  interval  (or  step  si/e)  and  that  the  current  step  value  is  A (i.e.  output  of  the  U.1  block 
is  1(A)).  The  following  configuration  statements  show  how  r may  be  represented  at  each  of  the 
j prcviinis  two  steps: 


.1 

V 5 

f(A-2) 

5 

U 6 

i-(A  i) 

6 

U 4 

y{k  1) 

4 

U 2 

,t(A-J) 

s 

U.1  9 

It  A ) 

6.2  Integrator 

Consider  the  first  order  differential  equation 

(ly  lit  =/  [r,  I'l 

Then,  for  the  second  order  Runge  Kutta  (i.e.  Modified  Euler)  method  used  in  the  integrator  (I) 
block,  the  dependent  variable  r.  which  is  the  output  value  of  the  I block,  is  defined  as 

r(f  h 2)  lit)  i (/i  2)/[r.,v(r)l 

r(r  /i)  = _f(f)  • A/(f  /i/2,  ,v(f  ■ /i  2)1 

where  li  is  the  time  step.  It  can  be  seen  that  at  each  half  time  step,  y is  defined  in  terms  of  quan- 
tities calculated  prev  iously.  It  is  therefore  evaluated  at  the  beginning  of  each  half  time  step  before 
execution  of  the  sort  list  (see  Section  .1.1).  Like  most  other  blocks  though,  an  I block  is  entered 
into  the  sort  list  only  when  its  inputs  have  already  been  sorted.  On  executing  an  I block  in  the 
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sort  list,  the  runclion  / IS  taleulaled  so  that  r may  be  derived  at  the  bcftinmng  of  the  next  ha'f 


time  step  ;is  described  above. 

The  first  order  l;ig  (11 ) block,  described  in  Section  6..1,  uses  the  same  integration  procedure 
as  the  I block  I hese  blocks  together  are  numbered  in  the  order  in  which  they  appear  m the  sort 
><  list,  this  number  is  placed  m the  integer  parameter  M for  each  of  these  bhx'ks.  thus  providing 

an  order  for  the  evaluation  of  r prior  to  entry  to  the  sort  list  at  each  half  time  step. 

I he  facility  to  reset  i to  its  initial  condition  or  hold  it  at  its  current  value  has  been  included 
in  the  I bloek  by  using  the  values  X2  and  X.f  as  switches’  (sec  Table  1).  Tor  the  reset  operation, 
the  value  of  i calculated  prior  to  execution  of  the  sort  list  is  overwritten  with  the  reset  value. 

It  may  therefore  be  important  to  ensure  that  the  integration  block  will  be  sorted  before  being 
used  as  an  input  to  another  block  (see  Section  .'.2).  Tor  the  hold  operation,  the  value  of  ciy  dt 
Is  set  equal  to  zero  so  that  the  current  value  of  i is  unaltered. 

(>.?•  First  Order  Tag 

I he  output  of  the  first  order  l;ig  ( Tl  I block  is  the  solution  of  the  dilferential  equation  ' 

I P2  d\di  XI  X2  X.'l 

for  ,in  initial  value  of  i equ,il  to  I’l.  Txcept  for  the  hold  and  reset  capability,  which  the  Tl  block 
does  not  have,  1 I blocks  are  treated  in  exactly  the  same  way  as  I blocks;  they  are  sorted  and 
executed  in  the  same  way  and  use  the  same  method  of  numerical  integration.  The  use  of  the  asso- 
ciated integer  parameter  \1  is  described  in  Section  6 2. 

7.  TSFK-DFFINFI)  HI  <)(  k.S 

C omplex  algebraic  expressions,  which  may  include  special  functions  such  as  trigonometric 
or  exponential  functions,  can  only  be  practicably  handled  in  a ’block  oriented’  language  by 
use  I’f  iiser-defmed  blocks  In  ( S\1P-I0(  \RI  ).  these  blocks  have  therefore  been  introduced 
;ind  ,ire  ideniilied  in  the  configuration  statements  by  a block  type  of  the  form  Um,  where  m is 
,in  integer  not  exceeding  fifteen;  this  number  may  be  increased  by  simple  coding  alterations  to 
the  mi'delling  program  I ach  block  may  have  up  to  three  parameters  and  three  sorted  block 
inputs  (as  deefired  m the  configuratii’n  statement),  but  may  have  any  number  of  unsorted  block 
inputs  When  reference  is  made  to  a block  of  type  Tm,  the  quantities  required  to  be  calculated 
must  be  delined  m a I OR  I R \N  subroutine  named  L^SI  Rm.  However,  each  of  these  sub- 
routines may  refer  to  other  (OR  I RAN  subroutines  or  functions.  .Ml  ’user-defined'  sub- 
routines and  ’other’  associated  subroutines  and  functions  are  compiled  and  loaded  into  core 
with  the  compiled  mi'delling  program  (see  Section  4 11). 

,\s  well  .IS  the  output  value  X of  the  user-defined  block  (referred  to  as  the  ‘principal’ 
output),  addition.il'  outputs  may  be  specified  within  the  subroutines  by  use  of  ’user  output’ 
blocks.  I hese  are  identified  in  the  conliguration  statement  by  a block  type  ‘L'O’  and  an  input 
HI  that  IS  the  associ.ited  iiser-delined  block  number.  They  are  not  entered  in  the  sort  list,  but 
are  listed  following  their  associated  user-defined  block  when  storing  the  model  (sec  Section  .^.2).  p 

Because  use  m.iy  he  mtide  in  these  subroutines  of  variables  stored  by  the  modelling  program 
in  labelled  COMMON  storage  areas,  a description  I'f  these  variables  in  given  in  Appendix  F. 

I I'r  an  understanding  of  the  rest  of  this  section,  the  reader  should  refer  to  the  subroutine 
example  given  in  Appendix  Ci  A model  incorporating  this  subroutine  is  described  at  the 
end  of  this  section. 

I he  same  user-defined  block  iy(ve  may  be  used  in  more  than  one  configuration  statement. 

Within  user-defined  subroutines  therefore,  it  is  advisable  to  avoid  referring  to  specific  block 

numbers  for  both  the  principal  and  additional  outputs.  This  may  be  achieved  in  the  following 

m. inner.  W hen  reference  is  made  to  a iiser-delined  subroutine,  the  block  number  B of  the 

.ippropriate  user-vleline<l  block  is  siorerl  in  ’I.’,  which  is  one  of  the  subroutine  arguments  (sec 

Appendix  Ci)  Outputs  of  each  block  used  in  the  model  are  stored  in  the  array  ‘C  ” so  that  the 

principal  output  should  be  set  in  the  array  clement  C (1  ).  Additional  elements  should  then  be  I 

set  in  eonsccutivc  elements  ( (I  I ).  ( ( I • 2).  . . . C(l  NUO),  where  NL'O  is  the  number  of  I 
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these  additional  outputs.  .An  examples  of  a set  of  appropriate  configuration  statements  (for 
NL'O  .^)  is  as  follows: 


12 

U2  8 5 

Principal  output  of  subroutine  USER2 

13 

UO  12 

First  additional  output  of  subroutine  USER2 

14 

UO  12 

Second  additional  output 

15 

UO  12 

Third  additional  output 

L'ser-defmed  subroutines  enable  the  user  to  overwrite  the  output  X.  inputs  Bl,  B2,  B3,  and 
parameters  F’l.  P2.  of  any  block.  As  an  important  safeguard  m using  UO  blocks,  it  is  there- 
fore important  to  test  for  compatibilits  between  the  user-detined  subroutine  and  appropriate  L'O 
configuration  statement  so  that  the  output  of  another  block  is  not  inadvertently  overwritten. 

for  UO  block  number  N.  the  tests  are: 

( 1 ) the  block  is  defined  as  a UO  block,  i.e.  Ml  R\(5*N  4)  } (the  type  number  of  a UO 

block,  see  Table  I ):  and 

(2)  Bl  is  the  associated  user-detined  subroutine  block  number,  i.e.  MTRX(5*N  — 3)=  L 
(i.e.  block  number  B). 

These  tests  need  only  be  performed  at  the  initial  time  step.  i.e.  when  TEST(5)  = I (see 
Appendix  F).  and  are  performed  in  subroutine  UOTFST  (see  Section  9. 1 ).  which  is  part  of  the 
basic  modelling  program.  For  NUO  consecutive  UO  blocks,  the  appropriate  statements  are  given 
in  the  example  of  Appendix  G;  for  a non-consecutive  UO  block  number  N.the  test  statement  is 

IF(TFST(5).EQ.I)CAL1.  UOTEST(N.L.MTRX) 

The  first  four  FORTRAN  siatemenis  in  Appendix  G would  generally  be  included  in  all 
user-detined  subroutines,  with  the  appropriate  value  for  m in  the  subroutine  name  'USERm'. 
For  convenience.  Bl.  B2,  B3.  XI.  X2.  X3.  PI.  P2.  and  I'3  arc  obtained  in  .Appendix  G from  the 
appropriate  arrays  for  the  particular  user-defined  block. 

Though  only  three  inputs  and  three  parameters  are  allowed  in  each  user-defined  block,  for 
all  other  blocks,  the  output  X.  inputs  Bl.  B2.  B3.  and  parameters  PI.  P2.  P3  may  be  obtained 
from  the  arrays ‘C‘.  ‘MTRX',  and  ‘PAR'  (see  Appendix  G for  definition  of  arrays).  FTowever. 
a knowledge  of  the  sorting  process  (sec  Section  3.2)  is  necessary  to  determine  whether  the  block 
output  values  are  at  the  present  or  previous  halftime  step.  For  each  UO  block,  the  parameters 
PI.  P2,  P3  and  the  inputs  B2  and  B3  may  be  used,  but  the  latter  inputs  are  not  required  to  be 
sorted,  since  UO  blocks  are  not  entered  in  the  sort  list. 

W hen  large  user-detined  subroutines  arc  envisaged,  it  is  advisable  to  construct  a small  user- 
detined  subroutine  that  refers  to  other  normal  I ORI  R.AN  subroutines  and  functions.  Any 
block  number  changes  may  then  be  confined  to  the  small  user-defined  subroutine,  thus  avoiding 
time  consuming  recompilation  of  large  subroutines. 

The  user-detined  subroutine  in  Appendix  G enables  the  divider,  constant,  and  weighted 
summer  blocks  used  in  the  model  of  the  non-linear  spring  problem  to  be  replaced  by  a user- 
detined  block  of  type  Ul  (see  Figure  3).  In  addition,  a user  output  block  provides  an  additional 
output  whose  value  is  the  same  as  the  replaced  weighted  summer  block.  It  is  not  suggested  that 
the  above  replacement  is  necessarily  desirable;  it  is  made  solely  for  the  purpose  of  presenting 
a simple  example  of  the  use  of  user-detined  subroutines.  Given  that  the  model  for  the  non- 
linear spring  problem  has  been  set  up  and  that  the  integration  and  output  control  parameters 
have  been  specified  (i.e.  the  model  is  at  the  completion  of  Stage  6 in  Section  4.1),  the  appropriate 
modifications  are  given  below: 

*CO\ 

COSHCL  RATIOSS  : 
lil.K.  TYPE.  Bl.  B2.  B3 
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4,  LI.  48.  10;  V DBLE  DOT 

PREV  Bl.K  4 DELETED 

5.10,4  S B*(Y  DOT)-F(V) 

6 

PRE\  BEK  6 DELETED 

17 

PREV  BLK  17  DELETED 

♦PAR 

PARAMETERS  : 

BLK.  PI.  P2.  PJ 

4,  5.  2 

*FLN 

EL'.KCTIOSS  : 

BLK  AO. 

MODEL  COMPLETE 

♦ 

At  this  st.ige,  the  model  may  be  run  using  the  'GOE'  command;  the  results  obtained  will 
be  identical  to  those  obtained  previously. 

8.  SPECIAL  FACILITIES 

The  two  special  facilities  described  here  are  the  solution  of  an  implicit  equation  and  the 
representation  of  discrete  systems  such  as  difference  equations.  The  two  blocks,  wye  (Y)  and 
vacuous  (V).  required  for  the  solution  of  an  implicit  equation  can  only  be  used  together  for  that 
purpose,  whereas  the  three  blocks,  unit  delay  (U).  zero  order  hold  (Z).  and  time  pulse  generator 
(T).  required  in  discrete  system  simulation  may  be  used  separately  in  other  applications  (see 
Table  I for  their  definitions). 

8.1  Solution  of  an  Implicit  Equation 

An  implicit  equation  (e  g.  f — /(  i ).  where  / is  a function  of  f and  possibly  other  variables, 
e.g.  time,  that  arc  independent  of  t ) cannot  be  solved  directly  in  CSMP-IO(ARL)  since  the  output 
value  of  the  block  that  calculates  / is  required  to  be  its  own  input  signal;  the  resulting  loop 
would  thus  lead  to  a failure  of  the  sort  algorithm.  However,  by  using  Y and  V blocks  (see  Table  I), 
an  implicit  equation  may  be  solved.  The  V block  ‘breaks’  the  sort  loop,  and  the  Y block  performs 
the  numerical  iteration  process  for  calculating  r;  the  final  value  for  r is  stored  as  the  output 
value  of  the  Y block. 

Assuming  that  the  function  /is  calculated  in  a user-defined  block,  U6  say.  Figure  4 shows 
the  block  structure  necessary  to  solve  the  implicit  equation  .r=/(t).  The  appropriate  block 
numbers  shown  correspond  to  the  following  set  of  configuration  statements; 

.f  Y 2 7 .V 

2 U6  7 /(y) 

7 V 

Any  other  blocks  that  arc  required  to  calculate  / must  also  have  the  output  value  of  the  V 
bliKk  as  an  input  signal.  The  iterative  method'®  used  to  solve  the  implicit  equation  y =f(y) 
IS  described  in  Appendix  H.  Included  in  Appendix  H is  a table  of  all  the  quantities  printed  for 
each  iterative  cycle  when  using  the  debugging  facility  (see  Section  5.5..f). 
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It  should  he  realised  that  the  numerical  method  used  to  solve  an  implicit  equation  will  not 
achieve  convergence  for  all  implicit  equations.  Very  few  numerical  techniques  are  capable  of 
solving  all  problems  of  a particular  type.  The  method  used  in  CSMP-IO(ARL),  which  is  a sig- 
nificant improvement  on  the  somewhat  crude  method  used  in  CSMP-IO,  should  be  viewed  as  a 
useful  technique  that  will  m most  cases  allow  the  user  to  solve  an  implicit  equation  without 
having  to  program  his  own  numerical  method  in  a user-defined  subroutine. 

S.2  Discrete  System  Simulation 

Discrete  systems,  such  as  dilVerence  equations,  may  be  modelled  in  CSMP-IO(ARL)  by 
using  V,  Z,  and  U blocks  (sec  Table  I).  Consider  the  second  order  difference  equation 

r(/)t)  /'[  r(/n  I ),  r(/n  2)]  for /n  2,3,.,, 

where  the  values  r(0)  and  r(  I ) are  known,  and  the  fimction/is  specified  in  a user-defined  block, 
1'3  sav  CSMP-I()(,\R1  ) is  written  primarily  for  continuous  systems  with  time,  r,  as  the  indepen- 
dent vantible.  When  dealing  with  purely  discrete  systems,  the  time  step  may  be  regarded  as  a 
cycle  step;  this  is  achieved  by  setting  the  pulse  train  period  (in  the  T block)  equal  to  the  time  step, 
f or  simplicity  in  the  following  description,  these  quantities  will  be  assumed  to  be  one  second, 
which  corresponds  to  a unit  cycle  step. 

1 he  combination  of  a Z,  U,  and  Z block  in  scries  is  such  that  if  the  output  value  of  the  first 
Z.  block  IS  \ {m).  then  the  output  value  of  the  second  Z block  is  \iiu  — I ).  Similarly,  a Z,  U,  Z,  U, 
and  Z block  series  combination  may  represent  r(/n),  i(m  I),  and  i(/i;  — 2)  as  the  output  value 
of  the  first,  second,  and  third  Z blocks  respectively,  figure  5 shows  the  block  structure  necessary 
to  represent  the  above  difl'erence  equation.  The  appropriate  block  numbers  shown  correspond 
to  the  following  set  of  configuration  statements: 


4 

f 

1 

S PULSE  TRAIN 

Block  1 = t 

6 

Z 

10 

4 

.nni) 

7 

z 

12 

4 

y{m  1) 

S 

z 

13 

4 

y{m-2) 

12 

u 

6 

13 

Li 

7 

10 

U3 

S' 

7 

f[y(m-\).  1 

fhe  output  viilue  of  each  block  used,  in  execution  order,  is  described  in  Table  3 for  each 
half  cycle  step  up  to  the  second  complete  step,  following  the  initial  step  (i.e.  k = 0),  the  output 
value  of  each  U block  remains  the  same  as  the  Z block  that  follows  in  the  series  combination  of 


TABLE  3 

BLOCK  OLTPLT  VALUES  FOR  A SECOND  ORDER 
DIFFERENCE  EQUATION  EXAMPLE 


Block 

number 

Type 

Step  (A  ) 

0 

} 

1 

J 

2 

12 

U 

.1(1) 

1(2) 

r(2) 

,1(3) 

If  3) 

13 

LI 

,1(0) 

.i'(l) 

.1(1) 

,f(2) 

,i-(2) 

4 

T 

1 

0 

1 

0 

1 

7 

Z 

1(1) 

r(  1 ) 

1(2) 

,1(2) 

.1(3) 

S 

Z 

1(0) 

1(0) 

1(1) 

,r(  1 ) 

1(2) 

10 

U3 

1(2) 

.1(2) 

r(3) 

,r(3) 

,1(4) 

6 

Z. 

1(2) 

.r(2) 

1(3) 

.1(3) 

.1(4) 

32 


Z and  U blocks  (see  Table  3).  By  setting  the  U block  parameter  values  in  the  way  shown  in 
l igure  5 this  property  also  applies  at  the  initial  step.  In  order  to  obtain  the  correct  output  value 
of  the  first  Z block  (i  e.  block  number  6)  at  k 0,  r(2)  should  be  calculated  using  the  difference 
equation  and  used  as  the  parameter  PI  for  the  Z block. 

9.  I)KS(  RIPTION  OF  MODELLING  PROGRAM 

The  computer  program  structure  and  storage  requirement  are  briefly  described  below, 
f urther  information  may  be  obtained  by  referring  to  the  FORI  RAN  listing. 

9.1  Program  Structure 

file  computer  program  consists  of  a small  main  program  (MAIN.)  used  to  open  the  arrays 
that  are  later  expanded,  a number  of  major  subroutines,  various  serv'-e  subroutines  and  functions, 
and  a BLOCK  D.AF.A  subprogram  which  specifies  the  list  of  block  type  names  (see  Appendix  I). 

9.2  Storage  Requirement 

Because  arrays  are  expanded  to  the  required  size  by  the  modelling  program  (see  Section 
4.1  2).  It  may  be  necessary  to  know  the  storage  dependence  on  the  number  of  various  blocks 
specified.  W hen  loaded,  the  basic  modelling  program  (see  Section  4.1.2),  without  any  additional 
subroutines  specified  by  the  user,  requires  I2K  words  of  core  store  (IK  1024).  On  expanding 
the  arrays,  the  number  of  additional  words  required  is  given  approximately  by  the  formula 

13  NBLK  • 4 NINT  NDF.L.AY  ■ 31  NFUN 

where  the  FORTRAN  integer  quantities  NBLK.  NINT.  NDEl.AY.  and  NFTJN  arc  defined 

III  the  COMMON  statement  labelled  ‘SIZE’  in  Appendix  E and  are  the  quantities  typed  (in  the 
above  order)  by  the  user  in  Section  4.1.2, 

If).  ( ()N(  LI  dint;  remark.s 

A description  of  the  ‘block  oriented’  simulation  language  CSMP-10(ARL)  has  been 
given  The  language,  which  has  been  developed  from  CSMP-IO,  is  written  mainly  in  FORTRAN 

IV  for  a POP- 10  computer.  Three  major  improvements  have  been  made.  F'irstly,  ‘user-defined’ 
blocks  written  as  lORFR.AN  subroutines  have  been  included.  These  blocks  may  have  only 
three  sorted  inputs,  but  may  have  any  number  of  unsorted  inputs.  ‘Dummy'  or  ‘user-output' 
blocks  have  been  created,  whose  output  value  is  defined  in  the  appropriate  user-defined  sub- 
routine. thus  allowing  any  number  of  outputs  for  each  of  these  blocks.  Secondly,  arrays  necessary 
to  store  information  on  each  block  used  arc  automatically  expanded  to  the  size  determined  by 
the  user.  Thirdly,  the  language  has  been  divided  into  a modelling  program,  which  is  used  to 
perform  the  model  simulation  and  store  the  output  in  a binary  file  on  a specified  storage  device, 
and  an  output  program,  which  is  used  to  print  and  plot  the  character  conversion  of  the  binary 
file.  Although  the  modelling  program  is  designed  to  be  run  interactively  from  a Teletype,  it  may 
also  be  run  non-intcractively  during  Batch  or  On-Line  processing.  A large  numTicr  of  minor 
improvements  have  been  made,  many  of  them  being  corrections  to  errors  in  the  available  version 
of  C'SMP-10.  The  Teletype  messages  and  responses  have  been  completely  revised,  and  much 
greater  detail  h;is  been  provided  here  than  in  descriptions  of  previous  versions  of  the  langugae 
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APPENDIX  A 


Time  Pulse  Generator  Algorithm 

The  lirst  pulse  occurs  when  XI  > 0,  and  a pulse  train  of  period  PI  continues  until  XI  0 
(if  at  all),  but  is  restarted  if  subsequently  XI  > 0.  Consider  the  case  when  XI  ^ 0 at  r t,;  then 


X = 0ift 

~ I if  t|  • «PI  $ t - 
= 0 if  r,  • (n  f0-5)PI 


f,  ■ (n  + 0-5)Pl 
^ r < t,  ‘ (n  ■ 


)PI 


w here  ii  =0.  I , 


and 


_ - - provided  XI  does  not  become  negative.  The  parameter  P2  and  integer 

parameter  M are  used  in  the  time  pulse  generator  algorithm,  which  is  described  as  follows  (h  is 
the  time  interval  for  integration  defined  in  Section  3.1): 

(a)  Initial  entry 

M = 0;  go  to  ( I ) 

ib)  Subsequent  entries 

(1)  If  XI  ^ 0.  go  to  (2) 

X 0;  M = 0;  exit 

(2)  If  M = I . go  to  (3) 

X I ; M = I ; if  PI  ^ 0,  print  appropriate  diagnostic  message 
(see  Appendix  B)  and  exit;  P2  = (/i  — PI)  2 • OOOOl *min(/i,  PI);  exit 

(3)  If  P2  ■ 0.  go  to  (4) 

Change  X ( 0 or  I);  P2  = P2  0-5*PI ; go  to  (3) 

(4)  P2  P2  • exit 

The  quantity  min(/i,  PI)  is  equal  to  either  h or  PI.  whichever  is  the  smallest;  the  quantity 
0 (X)0l*min  (/i.  PI)  is  therefore  an  error  term  that  ensures,  when  using  floating  point  arithmetic 
in  a computer,  that  the  equality  for  t occurs  at  the  lower  limit  (see  definition  of  X above). 


APPENDIX  B 


i 

I 

i , 

I 


i 
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Modi'llint>  Program  Diagnostic  Krror  Messages 
lo  (.llstinguish  the  following  error  messages  from  the  PORI  RAN  compiler  diagnostics,  all 
the  messages  below  are  preceded  b\  an  exclamaiion  mark  ( !).  Integer  constants  are  represented  by 
n and  ‘s.e.'  is  used  to  abbre\i;iie  ‘self-explanatory'. 

' -I/  /7  \//’/'  K>  P/l  /l'>H  H)  /I'RO  IS  lil.K  n when  using  di\ ider  block 
' l/  /■/,A//’/'  ro  TAKi:  SQRT  Of  A 17:'  SO-  IS  Rl  f n —when  using  half  power  block 

' hi  K O R CASSOT  (lO  fO  11  Y since  in  bmar_\  form 

' HI  K I RfSfRYfl)  lOR  TIMf  s.e. 

' HI  R n II  tS  AS  II  I f(,  U.  HI.R  TYHf  ("type")  sec  Table  I for  legal  block  types 

' H!  K n IS  SOT  ( l O HI  R when  attempting  to  define  an  additional  output  in  a user- 

defined  subroutine 

' HI  R n IS  SOT  ( V / HI.R  s.e. 

' HI  R n SOI  t SfP  when  specifxing  output  blocks 

' COPT  TRIOR  TO  CAI.I.  TO  TXI’ASP  L S R TCOC, SIST P -re-load  the  program;  if  error 
persists,  keep  data  files  and  consult  originator 

' CO  \l  U iS P TRROR  invalid  program  command  (sec  Table  2u) 

.'  COM  \l  ISP  SO  I.OSifTR  \ U IP  for  command  used  in  earlier  version 

' COSIHiS  M RT.APy  RTAP  TROM  MOPTL  IT  TUT  s.e. 

'CO\ll(,S  SOI  COMTITTT  s.e. 

' TRROR  m RISC  MOPTl.  TSTCCIIOS  s.e. 

' TRROR  l\  HI.R  n invalid  configuration  statement  (sec  Table  I) 

.'  T \ TTSSIOSS  M RT.APY  SIT  extensions,  including  period,  should  not  be  included  in 
‘channel’  command 

' T tlAI  ATTTMTT  TO  RT  tP  SOS-T.\TSTTST  PSK  Tll.T  "filename". "extension" 
causes  program  to  crash:  may  occur  only  when  running  non-interactivciy 

' I IS  II  TIMT  MCST  SOT  HT  I.TSS  TH  IS  ISTTTU.  TI.MT  s.e. 

.'  II  \(  TS  AIRTIP)  RTIP  TROM  MOPTl.  IT  ////.'—s.e. 

' / ( SCTS  SOT  CO.MTI.TTT  s.e. 

' II I TO  M PTI  K T SA.MT  must  be  I.OGn,  where  n 1 to  1 2 (excluding  5).  or  TTY 


ISU'ROriiR  PARAM  SPIX'II  ICATIOX  s.e. 


' ISCORRPCP  RPCORD ••record"  *♦**  try  again 

/.V7'«,  V COMROI.  PAR  A. US  Mil:  DEI)  s.e. 

' l\Ti:U\  IXTERl  AI.  MI  ST  \()T  Hi:  GREATER  TIIAX  TIME  RAXGE  —s.e. 

/\  I (/.//)  Bl.K  XO.  s.e. 

.'  ITERATIOX  TIMTT  OE  n EXCEEnED  TOR  I <S  > BEKS  n & n AT  TIME  t - 
see  Appendix  H 

> IP  (HI  K n)  TO  HIR  n IS  XOT  PEllXED  s.e. 

IP  HI  KIR  10  HIT  n IS  \OT  I HE  ASSOCIATED  L SER-DEEIXED  BEK  n when 
ilie  ii'er  I'uipiii  conliguration  siatement  is  inicompatihle  with  the  corresponding  user-defined 
statetnent 

' IP  H:  tor  > HI  K n IS  XOT  A I BEK  see  figure  4 
' rOGUM  I MT  n (I  READY  /V  LSE  h>  another  channel 
' MAX  OE  n II  \CTS  EXCEEDI.D  w hen  specifs mg  function  coordinates 
U \ \ OI  n / HI  KS  EXCEEDED  w hen  >pecif\ mg  configuration  statements 
!/4.V  OE  n I A TJ  HI  AS  EXCEEDED  s.e 
•XnXOEnl  HI  AS  EXCEEDED  se 

' U l\  OI  iPIX  COORD  P URS  EXCEEDED  see  Section  3.1.3 

UODI  I.  SOI  ( OUPI  I TT  all  configuration,  parameter,  and  function  statements  must 
lirsl  he  set 

' UOKI.  BORA  TO  DO  attempt  tci  execute  mc'del  before  setting  integration  or  output  con- 
trol parameters 

' ,U7  CORRESPOX  Dl\(/  COXEKi  STATEMEXT  when  specifying  parameters 
' \0  ((>\in,S  l\  UODII  IP  EII.E  s.e. 

\0  II  ACTS  l\  MODE!  IP  HIE  s.e. 

' XO  PiRiUS  l\  UODII  IP  I HE  s.e. 

' XOT  I \<H  (ill  (ORE  TO  SET  I P iRRA  YS  when  expanding  arrays  in  •EXPAND' 

’ O P COMROI  P \R\  US  SEEDED  s.e. 

' PIRIUS  II  REIDY  REID  EROU  MODEI  IP  EII.E  -s.e. 

' PI  RIOD  (PI)  I OR  I HI  K n ML  ST  HE  GREATER  TIIAX  ZERO  -see  Appendix  A 

' PM  nun  ABSCISS)  'real  constant'  DOES  XOT  EXIST  when  attempting  to  delete 
a coordinate  pair 

3S 


SOKI  l All.l  Ki:  AT  in.K  11 


see  Seel  ion  .V2 


-V  y.S'l'A  \ I'KROK  typing  error  or  inappropriate  response  to  Teletype  message  typed  by  the 
program 

I Si:K-ni:i  l\t:n  srimornsh:  \0T  LOAPI:!)  IOR  BLK  n when  the  user  has  not 
loaded  the  appropriate  user-delined  subroutine  or  subroutine  used  for  user-defined 
interpolation  in  the  fiinetion  block 


APPENDIX  C 

L'scr-Defined  Interpolation  for  Function  Block 

The  modelling  program  allows  the  user  to  specify  up  to  three  separate  interpolation  or 
curse-fitting  (e  g.  using  least  squares  technique)  formulae  for  use  m the  function  (I-)  block.  The 
appropriate  formula  to  be  used  for  each  I block  is  determined  by  the  values  set  for  PI.  Corres- 
poniling  to  PI  —I.  2,  and  2.  the  formula  used  must  be  programmed  in  the  FORTRAN  [] 

subroutines  INI  RPI,  INTRP2,  and  INTRP.'i  respectively.  The  first  two  statements  of  these  | 

subroutines  (e  g.  for  IN  FRP2)  take  the  form  1, 

SlJBROUTINi:  INTRP2  (I  . C,  MTRX,  PAR.  > OUT.  XIN.  Ml  IRST,  Ml  AST.  F)  I 

OIMFNSION  ('(I).  MTRX(I),  PAR(l),  1(1)  I 

where  I.,  M FRX,  and  P.AR  are  defined  in  Appendix  G,  I 

YOU  F is  the  function  output  value  set  by  the  user,  ] 

XIN  is  the  independent  variable  (abscissa)  for  which  YOUT  is  required,  and  ; 

F is  an  array  containing  the  coordinate  pairs  .X(K).  ^'(K)  < 

such  that  X(K)  F(2*K  l)andY(K)  F(2*K)  ' 

for  K MFIRST.  M FIRST  < I Ml  AST 

For  each  F block,  the  parameter  P.)  is  used  to  store  MFIRST,  and  the  integer  parameter  M is 
used  to  store  the  number  of  coordinate  parrs,  i.e  Ml  AST  MFIRST  > 


t 


AIM*KM)1\  I) 

Block  Output  for  Non-I. incar  Sprini*  Problem 

(I)  r.ibiilar  Oulpul 


time 

BLK*  9 

BLK*  4 

BLK*  40 

BLK*  IB 

0.00C>0l 

-1 .0BE*0l 

2.00E*0I 

0.00E-0I 

-I.BBE»B2 

1 .00E»00 

-3.27E*00 

- 1 . 42E*00 

9. 36E*00 

• 1 . 1 6E*0I 

2.00E*00 

4. 50E»O0 

-6.4IE*00 

S. S2E*00 

2.  I0E*0I 

3>00E*00 

6. 06E*00 

-6.39E«00 

•2. 44E400 

3.  68E»ai 

4.00E*00 

1 .86E*00 

1 . I0E*00 

-4.62E*00 

3.7|E*BB 

S.00E*00 

-2.02E«00 

2.0IE«00 

-2.98E*00 

-4. I0E*00 

6>00E*00 

-3.68E400 

2.90E*00 

•2.03C-0I 

• 1 . 41 E«0l 

7.00E*00 

-2. 69E*00 

8.9SE-0I 

1 .84E*00 

•B. 1 SE«00 

! 8<00E*00 

-7.47E-0I 

-4.38E-0I 

1 .84E*00 

- I . 49 E* 00 

* 9.00E*00 

8.00E-0I 

-7.97E-0I 

1 . I9E*00 

|.60E»00 

! I.00E»0I 

1 .S8E*00 

-7.B2E-0I 

3. 74C-0I 

3. I6E*00 

t |.|0E*0l 

1 . 60E400 

-5.83E-0I 

-2.92E-0I 

3.20E*00 

I.20E*01 

1. I0E*00 

-I.83E-0I 

•6. 46E-0I 

2.20E«00 

' |.30E*0I 

4.  1 7E-0I 

I .04E*0I 

-6. 76E-0I 

8. 33E-0I 

l.40E*0l 

- 1 . 75E-0I 

2. 63E-0I 

•4.8IE-0I 

-3. 51E-0I 

1 .S0E*0l 

-5. 1 6E-01 

2.85E-0I 

- 1 .97E-0I 

- I .03E*00 

1 .60E«0I 

- 8.80E'0I 

2.09E-0I 

5. 66E-02 

- 1 . I6E*00 

1.70E*0I 

- 4. 37E-0 1 

9. 1 6E-02 

2.08E-01 

-8.74E-01 

I.R0E*0l 

-2.02E-01 

- 1 . 64E-02 

2.43E-0I 

-4.05E-0I 

1 .90E*0I 

1 . 96E'02 

-8. 35E-02 

I .89E-0I 

3.92E-02 

2.00E*01 

1 «fe2E-0l 

- 1 .02E-0I 

9.25E-02 

3.24E-0I 

(2)  fiiaphical  Oulpiit 


BLK#  9 I -|.00E*01  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  l.00C«ei 

BLK#  4 I -0.00E*0l  ••••••••••••*•••••••  2.00E«0I 


TIME  I 


0. 00E-01  ♦ 
l.00E*00  I 
2.00E*00  I 
3.00E«00  I 
4.00E»00  I 
5.00E*00  1 
6.00E»00  1 
7.00E*00  1 
B.00E»00  I 
9.00C»00  I 
l.00E*0l  I 
I.I0E«0I  I 
l.20E*0l  I 
l.30E*0l  I 
l.40E*0l  I 
I.S0E*0l  I 
l.60E«0l  I 
l.70E*0l  I 

1. B0E*ei  1 
l.90E*0l  1 
2.00E*0I  I 
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DcbucKin};  Output  for  Non-I.itu’ar  S|)rinj>  I’roblcm 


TIMF 

B 

T 

X 

X 1 

X2 

X3 

M 

PI 

P2 

P3 

2.PI000F-0I 

9 

I 

-9.6097E*00 

3. 7726E*00 

0.0000E-01 

0. 0000E-0 1 

1 

- I .0000E«0l 

0.0000E-0I 

0.0000E-0I 

2.0000E-0I 

10 

E 

-9.2975E+01 

-9.6097E*00 

0.0000E-0I 

0.0000E-0I 

1 1 

0.0000C-0I 

0.0000E-0I 

l.0000E*00 

2.0000E-0I 

1 7 

M 

-8. 5430E*0I 

3.7726E*00 

-9.2975E*01 

e.0000E-ei 

0 

2.0000E'»00 

1 .e000E*00 

0.0000E-ei 

2.0000E-0I 

A 

/ 

I . 7086E*0I 

-8. 5430E*ei 

-S.0e00E»00 

0.e000E-0i 

0 

0.0000C-0I 

0.0000E-0I 

0. 0000E-0I 

2.0000E-0I 

48 

I 

3.  7726C*00 

1 . 7086E*01 

0.0000E-0I 

0.0000E-0I 

2 

0.0000E-0I 

0.0000E-0I 

0.0000E-0I 

2. 5000E-0I 

9 

I 

-9.421 IE*00 

4.  6269E*00 

0.0000E-0I 

0.0000E-0I 

1 

- 1 .0000E«0l 

0.0000E-0I 

0.0000E-CI 

2. S000E-01 

1 0 

E 

-8.9580E40I 

-9. 42 1 1 E» 00 

0.0000E-0I 

0.0000E-0I 

1 1 

0.0000E-0I 

0.0000E-0I 

1 .0000E*00 

2. S000E-0I 

1 7 

U 

-8.0326E*0I 

4. 6269E*00 

-8.9580E*0I 

0.0000E-0I 

0 

2.0000E«00 

1 . 0000E*00 

0.0000E-0I 

2.5000E-0I 

4 

/ 

1 . 6065E*0I 

-8.0326E*0I 

- 5.0000E*00 

0.0000E-0I 

0 

0.0000E-0I 

0.0000F-0I 

0.0000E-0I 

2. 5000E-01 

48 

1 

4.6269E*00 

1 . 6065E»0I 

0.0000E-01 

0.0000E-0I 

2 

0.0000E-0I 

0.0000E-0I 

0.0000E-0I 

3. 0000E-0I 

9 

I 

-9. 1 470E+00 

S. 379 1 E*00 

0.0000E-0I 

0.0000E-0I 

1 

- 1 .0000E*0I 

0.0000E-0I 

0.0000E-0I 

3.0000E-0I 

10 

E 

-8. 4647E*0t 

-9. 1 470E*00 

0.0000E-01 

0.0000E-0I 

1 1 

0.0000E-0I 

0.0000E-0I 

1 .0000E»00 

3.0000E-0I 

1 7 

W 

- 7.  3888E*0I 

5. 379 1 E400 

-8. 4647E*0I 

0.0000E-0I 

0 

2.0000E*00 

1 . 0000E*00 

0.0000E-0I 

3.0000E-0I 

4 

/ 

1 . 4778E*01 

-7.3888E*01 

- b. 0000E*00 

0.0000E-0I 

0 

0.0000E-0I 

0.0000E-0I 

0.0000E-0I 

3.0000E-0I 

48 

I 

5. 3791 E*00 

1 . 4778E»0I 

0.0000E-0I 

0.0000E-0I 

2 

0.0000E-0I 

0.0000E-0I 

0.0000E-0I 
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APPKNOIX  F 

Modelling  Program  Common  Variables 


For  each  labelled  COMMON  statement,  the  1 OR  I RAN  variables  appear  in  the  order  in 
which  they  are  described;  e.g.  for  COMMON  storage  area  labelled  'OOEVIM',  the  COMMON 
statement  is  COMMON  ODFVIM  ONLJM,  ODI  V.  101  \'. 


COMMON 

label 

FORTRAN 

identifier 

Declared 

type 

Description 

BATCH 

BATCH 

logical 

BA  R H TRUE,  if  a file  named  ‘BOMMP. 

IN'  is  found  on  the  disk 
I .Al.SE.  otherwise 

BUFF 

NBUFF 

integer 

Last  word  number  filled  in  array  ‘BLIFFER’ 

BUFFER!  126) 

real 

Block  output  bulfer  which  is  written  intoa  binary 
file  only  when  NBUFF  ^ 127 

NFl  AG 

integer 

NIT  AG  1 prior  to  output  of  first  block  value 
0 otherwise 

XPC 

real 

Percentage  change  required  for  output,  divided 
by  UX)  (i.e.  relative  change) 

OF  BUG 

DEBUG 

logical 

DF'BUG  TRUE,  if  debugging  in  operation 

.FALSE,  otherwise 

OMAX 

real 

lipper  time  limit  for  debugging 

OMIN 

real 

Lower  time  limit  for  debugging 

1 OGINO 

integer 

LOGIND  O ' for  debugging  output  directly 

on  the  teletype 

FXIRA2 

TY(4S) 

real 

l.ist  of  block  type  names  (including  Ul  to  UI5) 

I II  F 

POFV 

integer 

Model  output  device  name 

IFIFNA(2) 

integer 

Model  input  filename 

Ol  11  NA(2) 

integer 

Block  output  filename 

PFII  NA(2) 

integer 

Model  output  filename 

FIRST 

FIRST 

logical 

FIRST  TRUE,  for  block  output  at  initial 

time 

FALSF,  subsequently 

I INFSP 

LINESP(I) 

integer 

IINI^SP(I)  ‘(IH)';  array  containing  format 

specification  to  correct  line  spacing  following 
the  change  in  the  FORTRAN  operating  system 
from  'FORSE'  to  'FOROTS' 

NF.WFIL 

NEWFII 

integer 

NEWFH  0 if  block  output  file  has  not  been 

written  onto 
- 1 otherwise 

NUMB 

NOO 

integer 

Number  of  U blocks  used 

NEO 

integer 

Number  of  I and  Tl  blocks  used 

NIG 

integer 

Number  of  F blocks  used 

NCON 

integer 

Number  of  K blocks  used 

Nl  1ST 

integer 

Total  number  of  blocks  used 

4.1 
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APPENDIX  F (continued) 


1 


( OMMON 
label 

FORTRAN 

identifier 

Declared 

type 

NOMBA 

NUMBA 

integer 

ODEVIM 

ONUM 

integer 

ODEV 

integer 

IDEV 

integer 

PDEVIM 

PNUM 

integer 

SIZE 

NBLK 

integer 

NINT 

integer 

NEUN 

integer 

NDELAY 

integer 

NPNTS 

integer 

TEST 

TESTI 

integer 

TEST2 

integer 

TEST3 

integer 

TEST4 

integer 

TESTS 

integer 

TEST6  integer 

TEST7  integer 

TESTS  integer 

TEST9  integer 

TIM  DT  real 

OTS2  real 

TTOT  real 

TZERO  real 

TITLE  TITLE(I2)  real 

TSAMP 
T 


Description 


Number  ol  blocks  to  be  output 

Block  output  U>gical  unit  number 
Block  output  device  name  i e ITXin 
Model  input  device  name 

Model  output  logical  unit  number 

Mavimum  hliKk  number 

Maximum  numK*r  of  I and  T I bliK'ks 

Maximum  number  ol  I bliKks 

Maximum  number  ol  1 bliK'ks 

Number  ol  coordinate  pairs  used  in  E bltx:ks 

TEST  I I il  model  input  tile  not  opened 

2 if  opened  but  not  completely  read 
l)  it  completelv  read 

Model  input  logical  unit  number  if  equal  to  5, 
Teletv  pe  used  as  model  input  channel 
Tl  ST2  2 it  functions  have  been  read  in 

1 I'thcrvMse 

Tl  S7  4 I il  no  I bliKks  in  the  model 

2 il  data  for  I bliK'ks  not  completely 
read  in 

' if  completely  read  in 
TEST*'  I for  initial  time  step 

2 for  each  half  time  step 

1 for  each  complete  time  step 
4 for  error  during  execution 

for  run  terminated  by  typing  "*" 
o for  run  terminated  by  a Q block 
Number  representing  IcKation  of  last  character 
in  input  record  read  in  when  EFST6  0 or 
TEST6  ’’2,  new  input  record  is  read  in 
TEST"'  I if  integration  parameters  not  read  in 

2 otherwise 

TESTS  I if  output  control  parameters  not 
read  in 
2 otherwise 

TESld  I if  sort  process  not  completed 

2 it  sort  priK'ess  successfully  com- 
pleted 

Integration  interval  h 
h 2 

Final  time 
Initial  time 

Title  (up  to  60  characters) 

Output  time  interval 
Time  value.  / 


V 


real 

real 
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Example  of  I ser- Defined  Subroutine 


SUBROUTINE  USERl  ( L, C. M TRX . PAR ) 
INTEGER  TEST.B1.B2.B3 
COMMON/ TEST/ TESTt  9 ) 

DIMENSION  C(  1 ).MTRX(  I )* PARC  1 ) 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


♦ LIS  THE  BLOCK  NUMBER  OF  THE  USER-DEFINED  BLOCKI  THE  SAME 

♦ subroutine  COULD  BE  USED  BY  OTHER  BLOCKS. 

♦ 

♦ C(N)  IS  THE  OUTPUT  OF  BLOCK  N AT  THE  CURRENT  TIMEl  C(L)  MUST 

♦ THEREFORE  BE  SET  EQUAL  TO  THE  PRINCIPAL  OUTPUT  OF  THE 

♦ USER-DEFINED  SUBROUTINE.  SINCE  TIME  IS  STORED  IN 

♦ BLOCK  1.  C(I)  IS  THE  TIME  VALUE. 


MTRX  IS  AN  ARRAY 
MTRX ( 5»N- A) 
MTRX( 5*N- 3) 
MTRXt  S*N-2) 
MTRXt  S*N- I ) 
MTRXt  S»N) 


SUCH  THAT  FOR  BLOCK  N 
= TYPE  NUMBER 
= B1 
= B2 
= B3 

= QUANTITY  DEPENDENT  ON  BLOCK  TYPE  AND 
INSERTED  BY  THE  MODELLING  PROGRAM 


♦ PAR  IS  AN  ARRAY  SUCH  THAT  FOR  BLOCK  N 

♦ PAR<3*N-2)  » PI 

♦ PAR(3»N-1)  = P2 

♦ PAR(3»N)  s P3 

♦ 

♦ TFSKS)  « 1 FOR  INITIAL  TIME  STEP 

♦ =2  FOR  EACH  HALF  TIME  STEP 

♦ =3  FOR  EACH  COMPLETE  TIME  STEP 

♦ =4  FOR  ERROR  DURING  EXECUTION 


♦ OTHER  ELEMENTS  IN  ARRAY  ’TEST'  ARE  NOT  REQUIRED  IN  THIS 

♦ EXAMPLE 

♦ 

NUO  = number  of  consecutive  UO  BLOCKS 
DATA  NUO/ I / 

TEST  FOR  COMPATIBILITY  OF  UO  BLOCKS  WITH  USER-DEFINED  BLOCK 
TEST(S)  IS  SET  EQUAL  TO  4 IN  SUBROUTINE  UOTEST  IF  INCOMPATIBLE 

I Ft < TESTt 5) . GT. I ) .OR. <NUO.lt. 1 ) )GO  TO  22 

DO  I 0 J=  1 ,NUO 

N»L*J 


10  CALL  UOTEST  (N.L.MTRX) 

I Ft  TESTt  5) . EO. 4) RETURN 


non  !»  o o o 


Al’PKNDIX  C,  (continui'd) 


INPUT  VALUES  AND  PARAMETERS  OP  USER-DEFINED  BLOCK 

BI»MTRX(5*L-3) 

B2»MTRX(5*L-2) 

B3»MTRX( 5*L- I ) 

XI. 0 
X2.0 
X3.0 

I F(81 .GT. 0>X I .CCBI ) 

IF(B2.GT.0>X2.C(B2) 

IF(B3.GT.0)X3.C(B3) 

PI. PARI 3.L-2) 

P2.PARI 3*L- I ) 

P3.PARI 3*L) 

OUTPUTS  OF  USER-DEFINED  AND  UO  BLOCKS 

A.  PI 

B. P2 
YDOT.XI 
F.X2 

C(L* 1 ).B*YDOT*F 
C(L).-C<L*I )/A 

RETURN 

END 
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AFPKNDIX  H 


ltrrali>e  Method  for  Solving  an  Implicit  Equation 

The  algorithm  used  to  solve  an  implicit  equation  i h i ) is  given  below  (refer  Figure  4 
and  definitions  of  V and  ^ blocks  in  Table  1),  where  n is  the  iteration  index  (i.e.  f(n)  is  the 
approximation  for  v after  n iterations  and  T(/i|  is  an  improved  value  of 
n 0 

T(/i)  i„  when  t t„ 

V at  I It  2 when  t i„ 

(1)  i(/i  I)  /(r('0) 

If  n 0,  go  to  (2) 

C(«  1)  \{n  1):  go  to  (?) 

(2)  u [ i'(n  ■ I)  |■('l)]  [ ri'i)  v('/  1)1- </  ("  I I 

f(/i  ■ 1 ) <iy(n)  ■ { I </•,'■("  ■ ' ) 

(?)lf  .v(n)  10™.  go  to  (4) 

If  fill  ^ 1 ) 10  ™ go  to  (.*i) 

(4)  If  fill  ■ I)  fin)  fill)  Fc,  go  to  (5) 

II  n • I 

If  II  > Ni.  print  appropriate  diagnostic  message  (see  Appendix  B)  and  go  to  (5) 

Go  to  ( I ) 

(5)  r )■(//  1 ),  exit 


All  the  quantities  calculated  above  are  stored  either  as  bloek  output  values  or  parameters 
of  the  blocks  used.  I hese  quantities  are  punted  at  each  iteration  when  using  the  debugging 
facility  (sec  Table  H 1 ). 


TABLE  H.l 


Description  of  V and  Y Blocks  C orresponding  to  Debugging  Output 


Block  Type 

M 

PI 

P2 

P3 

V 

fin) 

Nv 

fin-\) 

yin) 

Y 

fill  ■ 1)  \in  ■ 1)  T('i  • 1) 

n 

Ec 

Ni  or  a* 

<7 

* P2  = Ni  for  initial  and  final  iterations 


a otherwise 


APPENDIX  I 

(1)  Modelling  Program  Structure:  Major  Subroutines 


Namc(s) 

EITect 

Section(s)  for 
Reference 

MAIN 

Conlrols  the  How  of  the  program 

— 

CSMO 

Zeroes  appropriate  quantities 

— 

CSMI 

Reads  in  conliguration  statements 

3.1.1,  4.1.4 

CSM2,  CSM3 

Sorts  configuration  statements 

3.2 

CSM4 

Reads  in  parameter  statements 

3.1.2,  4.1.4 

rSM5 

Reads  in  function  statements 

3.1.3,  4.1.4 

C'SMh 

Outputs  model  specification  statements 

4.1.5,  5.1.1 

C SM7 

Reads  in  integration  parameters 

4.1.6 

C SMSA 

Reads  in  output  control  parameters 

4.1.6 

t SMIO,  CSMI  1 

Executes  the  model 

3.3,  4.1.7 

CSMX 

Outputs  block  values  into  binun  Hie 

4.2 

CSMI  3 

Enables  examination  of  output  at  final  time  value 

.5.5.2 

i 
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APPENDIX  I (continued) 

(2)  Modelling  Program  Structure:  Service  Subroutines  and  Functions 


Nanic(s) 

Ellco 

Section(s) 

for 

Reference 

C'HARSH* 

Used  b\  •PRNT^  P■  ami  'Bl  KIN'  to  shift  alpha- 
numeric characters  stored  in  a EORTRAN  variable 

CHGIM. 

When  model  input  is  completed  for  a particular  tile. 

used  to  revert  the  model  input  channel  to  the  Teletype 

5.1.2 

CPU  + 

Obtains  CPU  run  time 

— 

DING 

Sounds  bell  on  the  Teletype 

— 

l OF  CHK 

Checks  for  end  of  model  input  tile 

5.1.2 

I XPAND+ 

F.\pands  an  array  to  the  defined  number  of  elements 

4.1.2 

I (HI  C Kt 

Searches  for  a (ile  on  the  disk 

MNPL’T.  KINPUT 

Used  to  read  in  data  in  the  form  of  individual  characters 

— 

GTITI  li 

Reads  in  title 

4.1.4 

INBl  KS 

L'sed  with  ‘CSMSA'  to  read  in  output  block  numbers 

4.1.6 

INOUT.  GETNAM 

Processes  the  ’channel'  command 

5.1 

1.  ERROR 

Used  by  dummy  versions  of  user-defined  subroutines 
to  record  error  .ind  type  diagnostic  error  message 

4 1.1 

OPE  BE 

Completes  block  output  tile 

4.2 

OR( 

l\ed  to  determine  whether  .i  file  may  be  overwritten 

5.2.1 

our 

Outputs  debugging  information 

5.5.3 

PAC  K 

I’sed  to  pack  d.ita  into  a binary  tile 

4.2 

PRNTVP.  BE  KIN. 
B1  KOUT.  t IIKIR 

Used  with  ’CS\ir  and  ‘CSM6'  to  output  configura- 
tion statements 

— 

RETAIN 

1 or  blocks  requiring  initial  conditums.  sets  the  first 

parameter  equal  to  the  last  calculated  output  value  so 
that  evecution  may  be  continued  at  a later  stage 

5.5.5 

REX' 

Used  t('  determine  whether  a file  may  be  read  from  the 

disk 

5.2.2 

RIjN* 

Allows  the  User  to  run  another  program  (e  g.  TR  ANS) 

5.5.4 

Sr  T.  SE-TC  OM 

Outputs  titling  and  timing  information  and  block  labels 

4.2 

STORE.  E MIN 

Used  to  manipulate  function  data 

3 1.3,  4.1.4 

TAGSRC 

Searches  nunlel  input  file  for  a particular  type  of  state- 
ment 

— 

EIMOUT 

Outputs  ( PU  run  time 

— 

TT^  C EIK  + 

Used  to  check  if has  been  typed  during  a run 

4.1.7 

1 OTE  ST 

Tests  for  compatibility  of  UO  bliKk  with  user-defined 
block 

7 

ART  1 ibrarv  subrout 

Hies. 

+ Written  in  MACRO-K). 
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AERO  NOTE  362 
FIG  1 


(a)  Parameters  PI , P2,  P3  are  shown  left  to  right 
If  only  one  is  shown,  parameter  is  PI. 

If  two  are  shown,  parameters  are  PI  and  P2 

(b)  Output  values  XI , X2,  X3  of  blocks  B1,  B2,  B3  are  shown  top  to  bottom. 
If  only  one  is  shown,  it  is  XI. 

If  two  are  shown,  they  are  XI  and  X2 


BLOCK  DIAGRAM  NOTATION  (refer  Section  3.1) 


A 


d 


A 


d 


AERO  NOTE  362 
FIG.  2 


(a)  Block  diagram  (refer  Fig.  1) 


(b)  Differential  equation;  y = — [by  + f(y)] /a  (Integration  interval,  h = 0.1s) 

(c)  Parameters;  a = 5,  b = 2,  y^  = — 10,  Vq  ~ ® 

y _io  -8  -6  -4  -2  0 2 4 6 8 10 

(d)  Function  f(y).  _4  0 4 16  36  64  100 

(e)  Output:  y,  y,  y,  f(y)  from  t = 0 to  20s  (Interval  = Is) 


(Continued  on  next  page) 


AERO  NOTE  362 
FIG.  2 (cont.) 


(f)  Configuration  statements 


4 

/ 

17 

6 

y 

6 

K 

-a 

9 

1 

48 

y 

17 

W 

48 

10 

by  + f(y) 

48 

1 

4 

y 

10 

F 

9 

f(y) 

(g)  Parameter  statements 

9 -10 

17  2 1 

6 -5 


(h)  Function  statements 


10 

-10 

-100 

-8 

-64 

-6 

-36 

-4 

-16 

-2 

-4 

0 

0 

2 

4 

4 

16 

6 

36 

8 

64 

10 

100 

block  number 


coordinate  pairs 


NON-LINEAR  SPRING  PROBLEM 


BLOCK  DIAGRAM  FOR  A SECOND  ORDER  DIFFERENCE  EOUATION  EXAMPLE 
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